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Abstract 

We present an investigation of the pion's electromagnetic form factor F W (Q 2 ) in the spacelike 
region utilizing two new ingredients: (i) a double-humped, endpoint-suppressed pion distribution 
amplitude derived before via QCD sum rules with nonlocal condensates — found to comply at the 
la level with the CLEO data on the transition — and (ii) analytic perturbation theory at the 
level of parton amplitudes for hadronic reactions. The computation of F 7r (Q 2 ) within this approach 



is performed at NLO of QCD perturbation theory (standard and analytic), including the evolution 
of the pion distribution amplitude at the same order. We consider the NLO corrections to the 
form factor in the MS scheme with various renormalization scale settings and also in the ay- 
scheme. We find that using standard perturbation theory, the size of the NLO corrections is quite 
sensitive to the adopted renormalization scheme and scale setting. The main results of our analysis 
are the following: (i) Replacing the QCD coupling and its powers by their analytic images, both 
dependencies are diminished and the predictions for the pion form factor are quasi scheme and 
scale-setting independent, (ii) The magnitude of the factorized pion form factor, calculated with 
the aforementioned pion distribution amplitude, is only slightly larger than the result obtained with 
the asymptotic one in all considered schemes, (iii) Including the soft pion form factor via local 
duality and ensuring the Ward identity at Q 2 = 0, we present predictions that are in remarkably 
good agreement with the existing experimental data both in trend and magnitude. 
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I. INTRODUCTION 



It is the purpose of this paper to review and discuss questions relating to the calculation 
of the electromagnetic pion form factor with an improved pion distribution amplitude (DA), 
derived from QCD sum rules with nonlocal condensates and to use QCD Analytic Per- 
turbation Theory (APT) @, beyond the leading order (LO). Before going into the 
details of this framework, let us expose, in general terms, what these two ingredients mean 
for the analysis and also make some introductory remarks. 

Hadronic form factors are typical examples of hard-scattering processes within QCD 0, 
UlEillll 11, 13 and clearly the first level of knowledge necessary to understand the structure 



of intact hadrons in terms of quarks and gluons. Such processes have been much explored 
both because of their physical relevance, as being accessible to experiments, and because 
the y a llow to assess nonperturbative features of QCD (for reviews, see, for instance, [3, 13 
EiL In the following, the discussion is centered around the pion's electromagnetic 

form factor. At a more theoretical level, "hard" means that at least some part of the process 
amplitude, recast in terms of quarks collinear to hadrons (in an appropriate Lorentz frame), 
should become amenable to perturbation theory via factorization theorems on account of the 
hard- momentum scale of the process, say, Q 2 , that should suppress factorized infrared (IR) 
subprocesses, thus ensuring short-distance dominance. Under these circumstances one can 
safely evaluate logarithmic scaling violations by means of perturbative QCD (pQCD) and 
the renormalizat ion-group equation. When no hard momenta flow on the side of the initial 
(incoming) or the final (outgoing) hadron, factorization fails and a renormalization-group 
analysis cannot be made, so that in order to calculate the non-factorizable part of the pion 
form factor, one has to resort to phenomenological models (prime examples of which are 
[rii Eil . licl l2l|). or employ theoretical concepts like the (local) quark-hadron duality 22| 
and their descendants |23j, |2J, |25( . 

Despite dedicated efforts in the last two decades, exclusive processes have failed to deliver 
a complete quantitative understanding within QCD for a variety of reasons, among others: 

• Limited knowledge of higher-order perturbative and power-law behaved (e.g., higher- 
twist) corrections to the amplitudes. 

• Presence of singularities (of endpoint, mass, soft, collinear, or pinch origin) that may 
spoil factorization in some kinematic regions. 

• Insignificant knowledge of hadron distribution amplitudes owing to the lack of a reliable 
non-perturbative approach. 

• Non-factorizing contributions that are not calculable within pQCD and hence intro- 
duce a strong model dependence. 

While it may still be not possible to clarify all these theoretical issues conclusively, we 
believe that significant progress has quite recently been achieved in understanding the pion 
structure both from the theoretical side — perturbatively j26[ 27, 28l I29I l30t and nonpertur- 
batively [H, HI, 3^, 3j| 34| — as well as from the experimental side |35Ll36l.l37j and associated 
data-processing techniques 38[ 3^, 4(j , a progress that could brin g to a cleaner com par ison 
between data and various theoretical QCD predictions H HIJ B H Q MM, H 
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EiL iif. 1 Moreover, a program to compute the electromagnetic and transition form factors 
of mesons on the lattice has been launched by two collaborations, 51, 5^ and [H^, that may 
provide valuable insights when it is completed. This situation prompts an in-depth review 
and update of these issues, in an effort to consolidate previous calculations of the pion's 
electromagnetic form factor and narrow down theoretical uncertainties. 
We will focus our present discussion on two main issues: 

(i) How QCD perturbation theory can be safely used to make predictions in the low- 
momentum regime where conventional power series expansions in the QCD cou- 
pling break down and nonperturbative effects dominate. Such an extension is 
based on recent works on "analytization" of the running strong coupling a s (Q 2 ) 
|E IE IE HE HE HE S2> EH ( see a l so HE| for a slightly different approach) and their 
generalization to the partonic level of hadron amplitudes, like the electromagnetic and 
the pion-photon transition form factor, or the Drell-Yan process, beyond the level of a 
single scheme scale (28 , IE HE iDJ • I n contradiction to the usual assumption of singular 
growth of a s (Q 2 ) in the IR domain, the QCD coupling in this scheme has an IR-fixed 
point, with the unphysical Landau pole being completely absent. In conventional per- 
turbative approaches, a choice of the renormalization scale in the region of a few Aqcd, 
as required, for instance, by the Brodsky-Lepage-Mackenzie (BLM) scale-fixing pro- 
cedure jH^, would induce singularities thus prohibiting the perturbative calculation 
of hadronic observables. Using APT, these singularities are avoided by construction, 
i.e., without introducing ad hoc IR regulators, e.g., an effective gluon mass 
and therefore the validity of the perturbative expansion (in mass-independent renor- 
malization schemes) is not jeopardized by IR-renormalon power-law ambiguities. In 
addition, APT provides a better stability against higher-loop corrections and a weaker 
renormalization-scheme dependence than standard QCD perturbative expansion 
El and Sec. IVITh . 
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(ii) How to improve the nonperturbative input by employing a pion DA which incorporates 
the nonp erturbative features of the QCD vacuum in terms of a nonlocal quark conden- 



sate OES 

EE EE Elf. This accounts for the possibility that vacuum quarks can flow 
with a nonzero average virtuality A^, in an attempt to connect dynamic properties of 
the pion, like its electromagnetic form factor, directly with the QCD vacuum struc- 
ture (we refer to for further details). Within this scheme, the pion DA (termed 
BMS [lj in the following) turns out to be double-humped with strongly suppressed 
endpoints (x = 0, 1), the latter feature bein g re lated to the nonlocality parameter X 2 ,. 
It has been advocated, for example, in |28|, |29| (see also jig, |69|), that a suppression 
of the endpoint region (which is essentially nonperturbative) as strong as possible is 
a prerequisite for the self-consistent application of QCD perturbation theory within a 
factorization scheme. 

In a recent series of papers plEEEEh two of us together with S. V. Mikhailov have con- 

O data [35j on the pion-photon transition using attributes from 
3 it |38l| . NLO Efremov-Radyushkin-Brodsky-Lepage (ERBL) 
Til . and detailed estimates of uncertainties owing to higher- 



ducted an analysis of the CL 
QCD light- cone sum rules 
U EE EE E2 evolution [re 



twist contributions and NNLO perturbative corrections |3Cj. These works confirmed the 



1 For theoretical predictions on F 7V (Q 2 ) in the timelike region, see, for instance 
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gross features of the previous Schmedding-Yakovlev (SY) analysis [38]; notably, both the 
Chernyak-Zhitnitsky (CZ) []~3j | pion DA as well as the asymptotic one are incompatible 
with the CLEO data [35[ at the 4cr- and 3<r-level, respectively, whereas the aforementioned 
BMS pion DA, which incorporates the vacuum nonlocality, is within the la error ellipse. 
Moreover, this approach revealed the possibility of using the CLEO experimental data to 
estimate the value of the QCD vacuum correlation length A" 1 . Indeed, it turns out that 
the extracted value A 2 ~ 0.4 GeV 2 is consistent with those obtained before using QCD sum 
rules (72! . FF3L [74| and also with numerical simulations on the lattice j7|| [76| . In addition, 
it was shown |4£| that the value of the inverse moment (x -1 ) w (/i 2 ) = L <fn(x; ii 2 )x~ l dx of 
the pion DA, calculated by means of an independent QCD sum rule, is compatible with 
that extracted from the CLEO data. These findings give us confidence to use the BMS 
pion DA (including also the range of its intrinsic theoretical uncertainties) in order to derive 
predictions for the electromagnetic pion form factor within the factorization scheme of QCD 
at NLO, presenting also results which include the non-factorizing soft contribution j^l |25| 
to compare with available experimental data. 

The structure of the paper is as follows. In Sec.|n]we shall recall the QCD factorization of 
the pion's electromagnetic form factor. Sec. IIHI deals with the basics of the pion distribution 
amplitude and its derivation from QCD sum rules with nonlocal condensates. The perturba- 
tive results for the pion form factor at NLO order, on the basis of the results given in [26j . are 
summarized in Sec. II VI whereas issues related to the setting of the renormalization scheme 
and scale are discussed in Sec. |V] The important topic of the non-power series expansion of 
the pion form factor in the context of Analytic Perturbation Theory is considered in Sec. I VII 
Our numerical analysis and the comparison of our results with available experimental data 
is presented in Sec. I VI II Finally, in Sec. IVIIII we give a summary of the results and draw 
our conclusions. Important technical details of the analysis are supplied in five appendixes. 



II. QCD FACTORIZATION APPLIED TO THE PION FORM FACTOR 

The outstanding virtue of factorization is that a hadronic process can be dissected in 
such a way as to isolate a partonic part accessible to pQCD. Provided that the partonic 
subprocesses are free of IR-singularities, then at large momentum transfer F n (Q 2 )/ f 2 ~ 
1/Q 2 , modulo logarithmic corrections due to renormalization. Hence, the amplitude for 
the electromagnetic pion form factor is short- distance dominated and can be expressed in 
terms of its constituent quarks collinear to the pion with the errors of this replacement being 
suppressed by powers of l/Q. Even more, one can rigorously dissect the QCD amplitude 
into a coefficient function, that contains the hard quark-gluon interactions, and two matrix 
elements corresponding to the initial and final pion states of the leading twist operator with 
the quantum numbers of the pion according to the Operator Product Expansion (OPE). In 
this way, one establishes that the coefficient function will scale asymptotically according to 
its dimensionality modulo anomalous dimensions controlled by the renormalization group 
equation. 

The pion's electromagnetic form factor is defined by the matrix element 

<7r+(P')l^(0)|7r+(P)> = (P + P'y^Q 2 ), (2-1) 

where J M is the electromagnetic current expressed in terms of quark fields, (P' — P) 2 = q 2 = 
—Q 2 is the photon virtuality, i.e., the large momentum transfer injected into the pion, and 
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is normalized to F^(0) = 1. Based on the above considerations, the pion form factor can 
be generically written in the form j^, 0, 11 , ^ 



F„(Q*) = FZ"*(Q 2 )+F, 



nnon-Fact 



(Q 2 



(2.2) 



where F^ act (Q 2 ) is the factorized part within pQCD and _p^ on - Fact (Q 2 ) is the non-factorizable 
part — usually being referred to as the "soft contribution" [23| — that contains subleading 
power-behaved (e.g., twist-4 and higher-twist) contributions originating from nonperturba- 
tive effects. It is important to understand that Eq. (|2.2j) becomes increasingly unreliable 
as Q 2 — > 0, owing to the breakdown of perturbation theory at such low momentum scales. 
Hence, we expect the real form factor to be different from the RHS of this equation at low 
Q 2 . We shall show in Sec. lVllDl how to remedy this problem. The leading- twist factorizable 
contribution can be expressed as a convolution in the form 



^Fact 



(Q 2 ; /4) = $l(x, (4) <8> T H (x, y, Q 2 ; /4, /4) ® $*(y, l4) 

x y 



(2.3) 



where ® denotes the usual convolution symbol {A{z)(&B{z) = L dzA(z)B(z)) over the 

z 

longitudinal momentum fraction variable x (y) and /xp represents the factorization scale 
at which the separation between the long- (small transverse momentum) and short- distance 
(large transverse momentum) dynamics takes place, with /% standing for the renormalization 
(coupling constant) scale. A graphic illustration of the factorized pion form factor in terms 
of Feynman diagrams is given in Fig. ^ 



• • ■ 
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FIG. 1: Illustration of the structure of the factorized pion form factor within pQCD at NLO of the hard 
scattering amplitude and the evolution effect of the pion DA. Hard gluons are indicated by broken lines, 
whereas the external off-shell photon is denoted by a cross. 



Here, Ts(x, y, Q 2 ; /x F , is the hard-scattering amplitude, describing short-distance in- 
teractions at the parton level, i.e., it is the amplitude for a collinear valence quark- ant iquark 
pair with total momentum P struck by a virtual photon with momentum q to end up again 
in a configuration of a parallel valence quark-antiquark pair with momentum P' = P + q 
and can be calculated perturbatively in the form of an expansion in the QCD coupling, the 
latter to be evaluated at the reference scale of renormalization /ij^: 

T K (x, y, Q 2 - /i 2 , /4) = «.(/£) if (x, y, Q 2 ) + (x, y, Q 2 ; /4) + . . . (2.4) 

The explicit results for the hard-scattering amplitude in LO and NLO accuracy are supplied 
in Appendix [XJ All transverse momenta below the factorization scale that would cause 
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divergences associated with the propagation of partons over long distances have been ab- 
sorbed into the pion DAs, which have the correct long-distance behavior, as dictated by 
nonperturbative QCD. 

Because the QCD perturbation series expansion in the strong coupling is only asymp- 
totic, this calculation bears an intrinsic error owing to its truncation that is of the order of 
the first neglected term ~ C/Q p , with C and p being, in general, dependent on the par- 
ticular observable in question — here the pion form factor. Lacking all-order results for the 
perturbative coefficients, one has to resort to fixed-order, renormalization and factorization 
scheme-dependent contributions to F^ act (Q 2 , fi 2 ^) that do not exceed beyond the NLO [26]. 
The truncation of this series expansion at any finite order introduces a residual dependence 
of the corresponding fixed-order or partly resummed hard-scattering amplitude and, conse- 
quently, also of the finite-order prediction for _Fj act , on the renormalization scheme adopted 
and on the renormalization scale chosen. In order that the perturbative prediction comes 
as close as possible to the physical form factor, measured in experiments — which is exactly 
independent of the renormalization (or any other unphysical scheme) scale — the best per- 
turbative expansion would be the one which minimizes the error owing to the disregarded 
higher-order corrections. This can be accomplished, for instance, by trading the conven- 
tional power-series perturbative expansion in favor of a non-power series expansion in terms 
of an analytic strong running coupling, performing the calculations in the framework of 
Analytic Perturbation Theory to be discussed in Sec. IVII Here it suffices to state that in 
this framework the QCD running coupling has an IR-fixed point and hence avoids eo ipso 
IR-renormalon ambiguities allowing to adopt a BLM scale setting procedure. 

By convoluting the finite-order result for the hard-scattering amplitude, expressed in the 
form of Eq. (|2.4|) . with the distribution amplitude (|3.5|) truncated at the same order in a s , 
an additional residual dependence on the factorization scheme and the factorization scale /xf 
appears. We show in App endix IB1 how to get rid of the factorization scale dependence at fixed 
order of perturbation theory (NLO) by proving that non-cancelling terms in F^ act (Q 2 ; fi 2 ^) 
are of order a^. For an alternative way of handling the /ip dependence, we refer the reader to 
For practical purposes, the preferable form of the convolution equation for F^ act (Q 2 ) 



is given by adopting the so-called "default" choice, i.e., setting in Eqs. (|2.3j) . (|2.4|) /ip = Q 2 . 
Note, however, that the same choice of scale in different schemes yields also to different 
results for finite-order approximants for the pion form factor 78] . Problems connected with 
heavy-quark mass thresholds in the (3 function are given below particular attention both in 
the hard-scattering part and in the evolution part. 

Another crucial question is whether the factorizable pQCD contribution to the pion form 
factor is actually sufficient to describe the available experimental data, or if one has to take 
into account the soft part as well. It has been advocated in 

B H S H H m M, 



that at momentum-transfer values probed experimentally so far, this latter contribution, 
though power suppressed because it behaves like 1/Q 4 for large Q 2 , dominates and mimics 
rather well the observed 1/Q 2 behavior. To account for this effect, we will include the soft 
contribution [25| (discussed in Sec. IV11 CJ) into our form-factor prediction when comparing 
with the data, albeit the poor quality of the latter at higher Q 2 makes it impossible to 
draw any definite conclusions about the transition from one regime to the other. Therefore, 
our main purpose in this paper is to calculate the factorizable contribution as accurately as 
possible. The calculational ingredients will be to 

• use as a nonperturbative input a set of pion DAs ^(x, /1q), derived in [l| from QCD 
sum rules with nonlocal condensates, with the optimum one, termed BMS model, 
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standing out; 

evolve (fi -Ax , n 2 ) by employing a kernel and corresponding anomalous dimensions up 
to NLO (7ol 71 , 80] both within the standard and the analytic perturbation theory; 



• em ploy a hard-scattering amplitude T H (x, y, Q 2 ; /z^) up to NLO order (2(1 8(J 81 
82I 83l . 84l 85^ . using both standard power and also non-power series expansions; 

• take into account the soft (non-factorizable) contribution, F^ on ~ Fa,ct (Q 2 ) , on the basis 
of the Local Duality (LD) approach when comparing the theoretical predictions with 
the experimental data. 

III. PION DISTRIBUTION AMPLITUDE 



A. Nonperturbative input 

Turning our attention now to the pion distribution amplitude, we note that ^ n (x, 
specifies in a process- and frame-independent way 2 the longitudinal-momentum xP distri- 
bution of the valence quark (and antiquark which carries a fraction x = 1 — x) in the pion 
with momentum P. At the twist-2 level it is defined by the following universal operator 
matrix element (see, e.g., [l3| for a review) 



(0 I d(z)^ l5 C(z,0)u(0) I vr(P)) 



z 2 =0 



iP" [ dxe lx(zP) ^ (x, fi 2 
Jo 



~z~ 2 ) 



/ $„.(x,^o) dx = /tt 
Jo 



with /„. = 130.7 ± 0.4 MeV |86] being the pion decay constant defined by 

(0|J(0) 7M 75«(0)|7r+(P)> = ipj n 

and where 

C{0,z) = Vexp 



-ig s t a Al(yW 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



is the Fock-Schwinger phase factor (coined the color "connector" in [87]), path-ordered 
along the straight line connecting the points and z, to preserve gauge invariance. The 
scale /2q, called the normalization scale of the pion DA, is related to the ultraviolet (UV) 
regularization of the quark-field operators on the light cone whose product becomes singular 
for z 2 = 0. 

Although the distribution amplitude is intrinsically a nonperturbative quantity, its evo- 
lution is governed by pQCD (a detailed discussion is relegated to Appendix EJ) and can be 
expressed in the form 



$ n {x, /4) = U (x, s; /4, ill) ® $ n (s, jil) 



(3.5) 



2 Provided the same factorization scheme is used for all considered processes 0, . 
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FIG. 2: Comparison of selectedpion DAs denoted by obvious acronyms: tp as (dotted line), <pcz (dashed 
line) 13], and </?bms (solid line) [J, defined by Eq. (|3.8() in conjunction with <|3.10|) . All DAs are normalized 
at the same scale /i§ « 1 GeV 2 . 



where Q n (s, /Iq) is a nonperturbative input determined at some low-energy normalization 
point /Xq (where the local operators in Eq. (j3.1j) are renormalized) — which is of the order 
of 1 GeV 2 — while U(x, s; /i 2 ,) is the operator to evolve that DA from the scale fi^ to the 
scale fip and is calculable in QCD perturbation theory. In the asymptotic limit, the shape 
of the pion DA is completely fixed by pQCD with the nonperturbative input being solely 
contained in /„.. 

Neglecting the dependence of the hard-scattering amplitude at large Q 2 , 3 it is conve- 
nient to introduce a dimensionless pion DA, ip n (x), normalized to 1, 

l4) = U <p*(x, nl) (3.6) 

that can be defined as the probability amplitude for finding two partons with longitudinal 
momentum fractions x and x "smeared" over transverse momenta < fi 2 , i.e., 

^{x^ 2 )= / [d 2 k ± ]i;(x,k ± ) , (3.7) 



o 



where [d kj_] is an appropriate integration measure over transverse momenta |12|, helicity 
labels have been suppressed, and a logarithmic pre-factor due to quark self-energy and 
photon-vertex corrections has been absorbed for the sake of simplicity into the definition of 
the pion wave function. 

The nonperturbative input, alias the pion DA at the initial normalization scale /i 2 ,, 
(p n (x, /ig), will be expressed as an expansion over Gegenbauer polynomials which form an 
eigenfunctions decomposition, having recourse to a convenient representation which sepa- 
rates the x and /i 2 dependence (a detailed exposition can be found in [3]). Then, the pion 
DA at the initial scale //q reads 

(p*(x, fil) = 6x(l - x) \l + a 2 (/j 2 ) Cl /2 (2x - 1) + o 4 (^) C 3 4 /2 (2x - 1) + . . .1 , (3.8) 



This actually means that for all initial kj_ i <^ Q 2 and analogously for all final 1^ ^ Q 2 , radiative 
corrections sense only single quark and gluon lines. 



s 



with all nonperturbative information being encapsulated in the coefficients a n . These co- 
efficients will be taken from a QCD sum-rules calculation employing nonlocal condensates 
0, HH , and we refer the reader to these works for more details. Here we only use the results 
obtained there. We found at /ig = 1.35 GeV 2 and for a quark virtuality of A 2 = 0.4 GeV 2 : 



a = 1 

a 6 = 5 x 10~ 3 



a 2 = 0.19 

a 8 = 4 x 10~ 3 



a 4 = -0.13 
aio = 4 x 10~ 3 



(3.9) 



One appreciates that all Gegenbauer coefficients with n > 4 are close to zero and can 
therefore be neglected. Hence, to model the pion DA, it is sufficient to keep only the first 
two coefficients, which we display below in comparison with those for the asymptotic DA 
and the CZ 13] model after 2-loop evolution to the reference scale /ig = 1 GeV 2 , i.e., 



Vbms 

^CZ 



i n = , n > 2 
i n = , n > 4 
i n = , n > 2 



a 2 

«2 



0.20 
0.56 



a 4 = -0.14 



/x 2 = 1 GeV 2 



(3.10) 



1 GeV 



The shapes of these DAs are displayed in Fig. |21 
At this point some important remarks and observations are in order. 

• The BMS pion DA, though doubly-peaked, has its endpoints x — > and x — > 1 strongly 
suppressed due to the nonlocality parameter A 2 Hence, fears frequently expressed in 
the literature that double-humped pion DAs should be avoided because they may 
emphasize the endpoint region, where the use of perturbation theory is unjustified, 
are unfounded. 

• The BMS pion DA approaches asymptotically <^ as in the endpoints from below, whereas 
V^cz approaches the asymptotic limit from above, which means that the endpoint 
behavior of the latter is dangerous until very large values of Q 2 . It is well-known 



that in the endpoint region 



x 



1 the spectator quark in the hard process, 
carrying the small longitudinal momentum fraction x, can "wait" for a long time until 
it exchanges a soft gluon with the struck quark to fit again into the final pion wave 
function. As a result, a strong Sudakov suppression |88| is needed in that case in order 
to justify the use of perturbation theory. In contrast, the endpoint behavior of the 
BMS DA is not controversial because, though doubly peaked, it does not emphasize 




FIG. 3: Percentage distribution (see text) of the first inverse moment in x of the BMS model DA 
comparison with the CZ one (a) and the asymptotic DA (b). 
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the endpoint regions. Even more, as Fig. |3] shows by plotting the first inverse mo- 
ment (x~~ 1 ) 7T , calculated as ' if 7r (x)x~ 1 dx and normalized to 100% (y-axis), the 
BMS DA receives in this region even less contributions than the asymptotic DA, as 
we explained above. 

By the same token, the Sudakov suppression of the endpoint region of the BMS DA is 
less crucial compared to endpoint-concentrated DAs. The implementation of Sudakov 
corrections using the analytic factorization scheme was considered in technical detail 
in (2^1 for the case of the asymptotic pion DA. Such an analysis for the BMS DA is 
more involved and will be conducted in a future publication. 

The deep reason for the failure of the CZ DA was provided in jH^, |6^, 8t|. The 
condensate terms in the CZ sum rules are strongly peaked at the endpoints x — > 
and x — > 1, the reason being that the vacuum quark distribution in the longitudinal 
momentum fractions is approximated by a ^-function S(x) and its derivatives. For that 
reason, the condensate terms, i.e., the nonperturbative contributions to the sum rule, 
force the pion DA to be endpoint-concentrated, with the perturbative loop contribution 
proportional tox(l—x) being insufficient to compensate these two sharp peaks at x = 
and x = 1. Allowing for a smooth distribution in the longitudinal momentum for the 
vacuum quarks, i.e., using nonlocal condensates in the QCD sum rules (as done in 
the derivation of the BMS pion DA), the endpoint regions of the extracted DA are 
suppressed, despite the fact that its shape is doubly peaked. 



Q 2 F 7 *^JQ 2 ) [GeVU- 




Q 2 [GeV 2 



FIG. 4: Light-cone sum-rule predictions for Q 2 F 7 « 7 _» 7r (Q 2 ) in comparison with the CELLO (diamonds, 
[9(1 |) and the CLEO (triangles, 01) experimental data evaluated with the twist-4 parameter value S^ W _ A = 
0.19 GeV 2 0E3- 

The predictions shown correspond to the following pion DA models: ipcz (upper dashed 
line) 13], BMS- "bunch" (shaded strip) [J, two instanton-based models, viz., |9l| (dotted line) and [j^f 
(dashed-dotted line), and the asymptotic pion DA ip as (lower dashed line) |9|,lllJ. A recent transverse lattice 
result [9^ is very close to the dash-dotted line, but starts to be closer to the center of the BMS strip for 
Q 2 > 6 GeV 2 . 



The endpoint behavior of the pion DA is the route cause why the pion-photon tran- 
sition form factor — which in LO is purely electromagnetic — calculated with the CZ 
pion DA was found j43] to overshoot the CLEO data. More recently, the analysis of 
the CLEO data using lieht-cone sum rules H H H E} has excluded the CZ pion 
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DA at the 4<r confidence level, while the BMS DA was found to be inside the la error 
ellipse (for A 2 = 4 GeV 2 ), whereas even the asymptotic DA was also excluded by the 
CLEO data at the 3a level. These quoted findings are reflected in the behavior of 
the predictions for the pion-photon transition form factor displayed in Fig. 0J which is 
based on the corrected version of 0] (the displayed strip is therefore slightly different 
from that in |47j). 

To make these statements more transparent, let us define the DA profile deviation factor 

A^^^ = l + a 2 + a 4 + ... (3.11) 

which quantifies the deviation of a model DA from the asymptotic one and supply its value 
in Table Ufor several pion DAs suggested in the literature in comparison with the constraints 
from the experimental data and theoretical calculations. Reading this Table in conjunction 
with Fig. 0J one comes to the conclusion that the BMS "bunch" provides the best agreement 
with the CLEO and CELLO experimental constraints, being also in compliance with various 
theoretical constraints and lattice calculations. 

TABLE I: Estimates for the Gegenbauer coefficients and the DA profile deviation factor A^, up 
to polynomial order 4 for the asymptotic, the BMS, and the CZ DAs compared with constraints 
derived from light-cone sum rules (LCSR), QCD sum rules with nonlocal condensates (NL QCD 
SR) for the DA and the inverse moment a and by analyzing the CLEO data. Also shown are 

the corresponding entries for instanton-based models (ADT, PPRWG, PR) and those associated 
with a transverse-lattice calculation — labelled Lattice. All values displayed are normalized at the 
scale fj? = 1.35 GeV 2 , corresponding to the scale of NL QCD SRs. 



DA models/methods 


<J2 


CI4 


a 2 + a 4 


a 2 — 04 




As 














1 


BMS 


0.19 


-0.13 


0.06 


0.32 


1.06 


CZ 


0.52 





0.52 


0.52 


1.52 


PPRWG [21] 


0.042 


0.006 


0.05 


0.04 


1.05 


PR [134] 


0.09 


-0.02 


0.07 


0.10 


1.07 


ADT [94] 


0.05 


-0.04 


0.01 


0.09 


1.01 


Lattice [93] 


0.08 


0.02 


0.10 


0.06 


1.10 


LCSRs [22] 


[0.07, 0.37] 








1.22 + 0.15 


NL QCD SRs for DA fl] 


[0.13,0.25] a 


[-0.04,-0.21] a 


[+0.02, +0.09] 


[+0.18, +0.46] 


1.06 ± 0.04 


NL QCD SR for (x" 1 ) [1] 






[+0.00, +0.20] 




1.10 + 0.10 


CLEO la-limits [40J 


[0.15,0.43] a 


[-0.60, -0.04] a 


[-0.21, +0.15] 


[+0.21, +1.00] 


0.97 + 0.18 


CLEO 2<7-limits [40J 


[0.11,0.47] a 


[-0.71,+0.07] a 


[-0.31, +0.25] 


[+0.07, +1.14] 


0.97 + 0.28 


CLEO 3<7-limits [40] 


[0.07,0.51] a 


[-0.82,+0.19] a 


[-0.41, +0.35] 


[-0.07, +1.28] 


0.97 + 0.38 



"Note that the uncertainties on the Gegenbauer coefficients 02 and 04 are correlated. Here, the rect ang ular 
limits of the fiducial ellipses extracted from the NL QCD SRs [|J and from the CLEO data in [3^, Etof are 
shown. 
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B. Perturbative NLO evolution 



Let us now discuss how the pion DA evolves at NLO using first standard perturbation 
theory to be followed by analogous considerations within APT. The evolution of the dis- 
tribution amplitude (|3.8j) proceeds along the lines outlined in Appendix O Taking into 
account only the first two Gegenbauer coefficients and LO evolution, one obtains 



x 



D,LO 



(// 



1} 



D,LO 







V ) (^x 



(3.12) 



where ci2 ' LO (/^f) and 04 ,LO (/ip) are given by (jC18|) taking recourse to (jC14|) . and D denotes 
"diagonal", while ND below stands for "nondiagonal" . On the other hand, the solution of 



2 

1.75 
1 . 5 

1.25 
1 

. 75 
. 5 
0.25 









//'/ 
lit S 

!'/ / 




\\\ 


- b / 


\0 — ''■r/V 


\ *\ 
\ \ 



X 




FIG. 5: The left panel shows the effect of the LO diagonal part of the evolution equation, Eq. I|3.12fl . on the 
BMS DA. The convex solid line denotes the asymptotic profile of the pion DA, the other solid one stands 
for (^bmsM a t 1 GeV 2 , while the broken lines represent Vm/is( x ) at ^> ^0, and 100 GeV 2 (with the larger 
scale corresponding also to the larger value of the DA at the middle point). Right panel: Comparison of 
<^bms ( E q- GUI)— solid line) and <^o ( Eq (|XT5|| — dashed line) at 20 GeV 2 to illustrate the effect of NLO 
evolution. 



the NLO evolution equation takes the form 



<fF°M) = ^' NLO (^/4) + ^^r NLO (*,/4) 



An 



(3.13a) 



where 



cpy h °(x,(x 2 F )=6x(l-x) 



1+ ^ a D,NLO ( ^ )C 3/ 2(2x _ i: 



(3.13b) 



n=2,4 



and 

^ D ' NL °(^, /4) = 6x(l - x)J2'aT^°(^) CT(2x - 1) . (3.13c) 

n>2 

The coefficients a° ,NLO (/i F ) and a^ D ' NLO (/i|) are given in ()C19b|) and (jC19c|) . respectively, 
by employing ()C15|) . while denotes the sum over even indices only. Note that, although 
the input DA, (f^x, //q), was parameterized by only two Gegenbauer coefficients ai and 04, 
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higher harmonics also appear due to the nondiagonal nature of the NLO evolution. 4 The 
effect of the inclusion of the LO diagonal part of the evolution kernel is important, as one 
sees from the left part of Fig. which shows this effect for the BMS pion DA. 

On the other hand, from the right part of Fig. |5J we deduce that the NLO nondiago- 
nal evolution is rather small. We note that in the above computation the exact two-loop 
expression for a s (95j in the MS -scheme (Aq CD = 410 MeV, Nf = 3) was employed, cf. 
(I6.15|) . in which matching at the heavy-flavor thresholds M4 = 1.3 GeV and M5 = 4.3 GeV 
(with Mi = M2 = M3 = 0) has been employed 96]. A discussion of the relation of this 
exact solution to the usual approximation, promoted by the Particle Data Group 0], has 
recently been given in (see also Appendix C). 



IV. PION FORM FACTOR AT NLO: ANALYTIC RESULTS 



The NLO results for the hard-scattering amplitude Th are summarized in Appendix [XJ 
Setting in ()2.4|) /ip = Q 2 and taking into account the NLO evolution of the pion DA ip via 
(13~T3J) . we obtain from (Q 



where the LO term is given by 

^ LO (Q 2 ;/4) = « s (/4)^ LO (Q 2 ) 



Q 2 ^°(Q 2 



^f 2 



l+a 2 {Q)+a 



D,NLO 
4 



(Q 2 



(4.1) 

(4.2) 
(4.3) 



and the calligraphic designations denote quantities with their a s -dependence pulled out. 
In order to make a distinction between the contributions stemming from the diagonal and 
the nondiagonal parts of the NLO evolution equation of the pion DA, we express the NLO 
correction to the form factor in the form 



77-NLO / ' ( n2. 2 \ a s(/%) 

IT {Q 1 /%J = 



7T 



[^ D ' NLO (Q 2 ; /4) + ^ ND ' NLO (Q 2 ; AW = oo 



and write the diagonal contribution 

^ D ' NLO (Q 2 ;/4) = &o^ 1 ' /3) (Q 2 ;4) + ^ 1 ' F) (Q 2 ) + ^ 1 ' G) (Q 2 ) 

as a color decomposition (in correspondence with (|A3jl ) in terms of 



(4.4) 



(4.5) 



5 3 
3 + " 



(43/6)a°' NLO (Q 2 ) + (136/15)a 4 D ' NLO « 2 ) 



1 , D.NLO 

1 + a 2 



(Q 2 ) + a?' NLO (Q 2 ) 



In 



x 



D,NLO/^,2\ , D,NLO/ /n 2> 



/4 



(4.6a) 



Since a„ ' (/ip) decreases with n, for the purpose of numerical calculations, we use an approximate 



form of ^JjJd (x, /ip) in which we neglect a. 



ND,NLOf,,2 



(f4) for n > 100. 
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Q 2 ^ ¥ \Q 2 ) = 2*fl 



71 



— a 



D.NLO 



(Q 2 



497 161 D . 



NLO 



36 



24 



(Q 2 



D,NLO 



(Q 2 ) 



1123 9793 D . 



NLO 



450 



300 



(Q 2 ) 



1387 
~50~' 



D,NLO 



(Q 2 



(4.6b) 



and 

Q 2 ^ G \Q 2 ) 



-0.67 + a 



D.NLO 



(Q 2 



18.70 + 16.35 a 



D.NLO 



(Q 2 



+ a 



D,NLO 



(Q 2 ) [24.23 + 36.76 a°' NLO (g 2 ) + 20.26 a4' NLO (g 2 )]} , 



(4.6c) 



where the superscripts F and G refer to the color factors Cp and Cq = Cp — Ca/2, respec- 
tively. Note that for the matter of calculational convenience, we also display the sum of 
these two terms (cf. Eq. (|A7|) ): 



g2 ^(i,FG) (g2) = 2 vr/ 2 {-15.67 



D,NLO 



(Q 2 



- a A 



D,NLO 



21.52 - 6.22 a 



(Q 2 ) 7.37 - 37.40 a°' NLO (Q 2 ) - 33.61 a 



D,NLO 



{Q 2 



D.NLO 



(Q 2 



(4.7) 



On the other hand, the nondiagonal term reads 



g2 ^ND,NLO (g2 .^ )=47r/ 2 l + a D.NLO (g2)+G D,NLO (g2) £ , 



ND.NLO 



(Q 2 ) , (4-8) 



k=2 



where iVMax denotes the maximal number of Gegenbauer harmonics taken into account in 
order to achieve the desired accuracy. 

As it was shown in Ref. jl^j], the effects of the LO DA evolution are crucial. In order to 
investigate the importance of the NLO DA evolution, we compare the predictions obtained 
using the complete NLO results, given above, with those derived by employing only the LO 
DA evolution via 1)3.12)1 . The corresponding expressions in this latter case follow from those 
above by performing the replacements 



D,NLO 



D,LO 



a"—-* a"— and a^ D ' NLO -> , (4.9) 
so that the contribution J r ^ T>,NLO (Q 2 ] fj^) is absent. Introducing the notation 

•F-k = 3~-k I LO DA evolution > 

(4.10) 

we analyze the relative importance of the various contributions (LO, NLO, and Local Duality 
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(LD) part — see Sec. IVH Cj) by defining the following ratios 

^° m + k/tt) j; nd ' nlo m n m ^ 



R (Q 2 , iV Max ) 
R (Q 2 , N Mslx ) 

-Rmod {Q 2 ) 
-Rmod {Q 2 ) 



^° (Q 2 ) + (a s /n)^ D > NLO (Q 2 ; AW « oo 



(Q 2 ) + « s ^ LO (Q 2 ) + (« 2 /tt)^ 



ND.NLO 



(Q 2 ;iVi 



Max) 



ND,NLO 



Max 



OO 



^ LO (Q 2 ) . 

^° (Q 2 ) + (a s /vr)^ ND ' NL ° (Q 2 ; A Max « oo) ' 
F^ D (Q 2 ) + c^ LO (Q 2 ) 

(Q2) + ag jrLO (Q2) + (Q 2 /7r):F ND,NLO (Q2 . ^ ^ ^ 



(4.11) 

(4.12) 
(4.13) 
(4.14) 




R (Q 2 , A" Max ) 



20 



80 100 



"Max 



40 60 80 

^Max 



100 



FIG. 6: Left: The ratio R (Q 2 , A^ax) , defined in Eq. Q4.11|l . for three different values of Q 2 as a function 
of iVMax- The blue dotted curve corresponds to Q 2 = 2 GeV 2 , the green solid curve to Q 2 = 10 GeV 2 , and 
the red dashed one to Q 2 = 50 GeV 2 . Crosses of the same color represent the values of R (Q 2 , A^Max = 4), 
whereas triangles refer to i? mo d (Q 2 )i Eq. I|4.13[) . Right: The same designations hold as for the left side, but 
now for the ratio R (Q 2 , Ajyiax) given by Eq. (|4.12|) and correspondingly for -R mo d — Eq. i|4.14|) . 

These ratios are displayed graphically in Fig. El for the BMS DA (a 2 (/io) = 0.2, a^/i 2 ,) = 
—0.14) in the region A^Max = 4 — 100. We infer from this figure that, adopting in our 
calculations an accuracy on the order of 99.5%, we can safely neglect the non-diagonal part 
of the NLO evolution equation and use for the pion form-factor computations to follow the 
approximate expression (omitting the superscript Approx) 



F Fact-App r ox (g2 . ^ = a ^ ^LO (g2 . ^ + 



a 2 (/x 2 



^ D ' NLO (Q 2 ;/4)- (4-15) 



7T 



Actually, the difference between ^°' NLO (Q 2 ; /4) and J r °' NLO (Q 2 ; /i|) is of the order of 
a s (Q 2 ), so that it is safe to use everywhere only the LO evolution. We have verified in 
our numerical calculations that the difference is indeed less than 1%. 



V. SETTING THE RENORMALIZATION SCHEME AND SCALE 

The choice of the expansion parameter represents the major ambiguity in the interpreta- 
tion of the pQCD predictions because finite-order perturbative predictions depend unavoid- 
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ably on the renormalization scale and scheme choice. 5 If one could optimize the choice of 
the renormalization scale and scheme according to some sensible criteria, the size of the 
higher-order corrections, as well as the size of the expansion parameter, i.e., the QCD run- 
ning coupling, could serve as sensible indicators for the convergence of the perturbative 
expansion. In what follows, we shall consider several scheme and scale-setting options. 



A. MS scheme 

The results we have presented in the previous subsection were obtained in the MS renor- 
malization (and factorization) scheme. Let us discuss the choice of the renormalization scale 
/xr in some more detail. We see that in our NLO results, Eq. (|4.4)1 . this dependence is con- 
tained in the coupling constant a s (/j R ) as well as in the NLO correction J r i 1,/3 \ The latter 
correction is proportional to the 60 coefficient of the (3 function and is A/-dependent. Hence, 
a natural question arises: How can we determine the right value of Nf in the form-factor 
expression? 

We propose here to apply the following procedures, 
(i) The first one concerns the standard choice /j r = Q 2 and suggests to shift /i R at the 
heavy-quark threshold in order to ensure the continuity of the form factor according to 

F^(Q 2 ;/ i R = Q 2 )L, for Q 2 <Mf; 



F^(Q 2 ) 



Q 2 + < M 5 
F^(Q 2 ; n 2 R = Q 2 + 6 ii\ + 6&\ for Mf < Q 2 + ^\ . 

As a result, we have to fulfill the following matching conditions 

F^(M 2 ; M 2 4 )\ Nf=3 = F^(M 2 ; M\ + 5^)\ Nf=4 ; (5.1b) 
F^(M 2 - 5/4; M 2 )\ Nf=4 = F^\M 2 5 - 6^ M 2 + */iJ)L /=B - (5-lc) 

(ii) The second procedure addresses specifically the BLM scale setting /i R = /x BLM . In this 
case, the only problem is the small value of the BLM scale (see Table |H} due to the fact 
that the bo-term is completely absent and Ay-dependent terms do not arise. Therefore, we 
propose to implement the BLM scale setting only above some minimal scale: /i m i n - Below 
this scale, which is in the range of the typical meson scales and hence only the light-quark 
sector (Nj = 3) contributes, we fix /i R = /i^ in and set Nf = 3 using the ^^^(Q 2 ', /i ^ in )-term 
in the form provided by (J5.13|) — more explanations will be given shortly. 

The truncation of the perturbative series to a finite order introduces a residual dependence 
of the results on the scale /jr, while the inclusion of higher-order corrections decreases this 
dependence. Nonetheless, we are still left with an intrinsic theoretical ambiguity of the 
perturbative results. One can try to estimate the uncertainty entailed by this ambiguity 



5 Actually, to the order we are calculating these dependencies, they can be represented by a single parameter, 
say, the renormalization scale, because bo and b\ are renormalization-scheme invariant. 



16 



(see, for example, |26j) or choose the renormalization scale /xr on the basis of some physical 
arguments. 

The simplest and widely used choice for /ir is to identify it with the large external scale, 
i.e., to set 

/4 = Q 2 , (5.2) 

the justification for adopting this choice being mainly a pragmatic one. However, physical 
arguments suggest that a more appropriate scale should be lower. Namely, since each ex- 
ternal momentum entering an exclusive reaction is partitioned among many propagators of 
the underlying hard-scattering amplitude in the associated Feynman diagrams, the physical 
scales that control these processes are related to the average momentum flowing through the 
internal quark and gluon lines and are therefore inevitably softer than the overall momen- 
tum transfer. To treat this problem, several suggestions have been made in the literature. 
According to the so-called Fastest Apparent Convergence (FAC) procedure j9?i HI[ , the scale 
/xr is determined by the requirement that the NLO coefficient in the perturbative expansion 
of the physical quantity in question vanishes which here means 

^r°(Q 2 ;/4 = /4 AC ) = o. (5.3) 



On the other hand, following the Principle of Minimum Sensitivity (PMS) [99|, llOOl . llOlL 
102], one mimics locally the global independence of the all-order expansion by choosing the 



renormalization scale /xr to coincide with the stationary point of the truncated perturbative 
series. In our case, this reads 

4[W;,I) + ™;,I)]L | = 4ms =o. (5.4) 

In the Brodsky-Lepage-Mackenzie (BLM) procedure , all vacuum-polarization effects 
from the QCD /5-function (i.e., the effects of quark loops) are absorbed into the renormalized 
running coupling by resumming the large (boas) n terms giving rise to infrared renormalons. 
According to the BLM procedure, the renormalization scale best suited to a particular 
process at a given order of expansion can be, in practice, determined by demanding that the 
terms proportional to the /3-function should vani sh. This naturally connects to conformal 
field theory and we refer the interested reader to 103] for a recent review. The optimization 



of the renormalization scale and scheme setting in exclusive processes by employing the BLM 
scale fixing was elaborated in jiBf and in references cited therein. The renormalization scales 
in the BLM method are "physical" in the sense that they reflect the mean virtuality of the 
gluon propagators involved in the Feynman diagrams. According to the BLM procedure, 
the renormalization scale is determined by the condition 

^W;/4 = /4lm)=0- (5.5) 

For calculational convenience, we express /x R in terms of Q 2 : 

/4 = <Q 2 ) Q 2 (5.6) 

and proceed to calculate this quantity in the above-mentioned scale-setting schemes. Then, 
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TABLE II: Scales //pms> MFAC) ^blm> and /iy for the asymptotic, the BMS, and the CZ DAs. 



DA Q7/4AC 




Q 2 1 /^PMS Q 2 1 ^BLM 


<?7a£ 


Q 2 


As 18 




27 106 


20 


any 


BMS 16 - 20 




24 - 29 105 - 117 


20 - 22 


1 - 50 GeV 2 


CZ 146 - 62 




217 - 92 475 - 278 


90 - 52 


1 - 50 GeV 2 


the FAC procedure leads to 








a F Ac(Q 2 ) = exp 


5 3 + (43/6)a°' NLO (Q 2 ) + (136/15)a°' NLO (Q 2 ) 


o l + a 2 \Q)+ a i \Q) 








4 ^ 1 ' FG) (g 2 ) + ^ ND ' NLO (Q 2 )" 

bo ^°(Q 2 ) 




(5.7) 


which can be related to the PMS procedure via 










a PMS (Q 2 ) = e~ Cl/ VAc(Q 2 ) 




(5.8) 


with ci = &i/&o- This value corresponds to the stationary point (the maximum) of the NLO 
prediction for _pj act . 

On the other hand, for the BLM scale one obtains 






/^BLM — a BLM(Q 2 ) Q 2 , 




(5.9a) 


where 










ciblm(Q 2 ) = exp 




5 3 + (43/6)a 2 3 ' NLO (Q 2 ) + (136/15)a?' NLO (Q 2 )~ 


(5.9b) 




o l + a 2 (Q ) + a 4 (Q ) 



The values of the scales /ipms, /^fac, and /zblm for the asymptotic, the CZ, and the BMS 
DAs, defined in (13.10)1 . are listed in Table ITU One notices that the BLM scale is rather 
low for all considered DAs. This makes its applicability at experimentally accessible Q 2 
values rather questionable. But it is possible to improve this scale-setting procedure in the 
following way. 

First of all, let us rewrite the BLM prescription in the more suggestive form 

faq 2 ^1),^\Q 2 ^ 2 r) = \^°(Q 2 )^ 

(5.10) 

{fAQ 2 ; v 2 BLM(Q 2 ))^i m (Q 2 ; vIlm(Q 2 )) = 0} . 

It becomes evident that when the BLM scale yields a s values close to unity, perturbation 
theory breaks down. To avoid this happening, one can, of course, introduce ad hoc a cutoff 
for a s , operative, say, above 0.5 — 0.6, or one can "freeze" a s at low Q 2 scales to some finite 



/4 



./^blmW 2 



BLM 
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value by introducing, an effective gluo n mass |45j, |63[. Still another possibility is to use the 
analytic coupling as done in |28l. I29I] (see next section). 

In order to protect the BLM scale from intruding into the forbidden nonperturbative soft 
region, where perturbation theory becomes invalid, one can make use of a minimum scale, 
/i min , based on t he grounds of QCD factorization theorems and the OPE, as applied for 
instance in [Tol EM [Tol, EflJ and also in 0: 

fj? min > nl. (5.11) 

Here /Zq stands for a typical nonperturbative (hadronic) scale in the range [0.4—1.5] GeV 2 and 
corresponds roughly to the inverse distance at which the parton and hadron representations 
have to match each other. Note that the smaller /^ in is chosen, the deeper the endpoint 
region x — > 1 can be explored for smaller values of Q 2 . It is intuitively clear that the 
typical parton virtuality in the (hard) Feynman diagrams — let us call it fi 2 — should not 
become less than its counterpart in the pion bound state: fi 2 . Because the latter is linked 
to the scale /ip, the scale should be limited from below by this scale. Consequently, 
we assume that the following hierarchy of scales — partonic (i.e., perturbative) and hadronic 
(i.e., nonperturbative) — holds: 

Ag < /^o < l^q < /^-scheme ■ (5-12) 

Then, if /i| LM < //q, one obtains instead of Eq. (|5.1U|) . the IR protected version (termed in 
our analysis BLM prescription) 

[q 2 ,/4lm(q 2 )] [Q 2 ^Ilm(Q 2 )] = 

(5.13) 



BLM 



2 N 



4 W ™ ' V/4lm(Q 

This modification of the BLM scale setting enables us to treat the problem of the Nf- 
dependence of the (3 function in the term !F^'^\Q 2 ] /j^) without any further assumptions or 
modifications. Because of the fact that the scale [i\ is now bounded from below by (J5.11j) . 
one is not faced with ambiguities related to the variation of the number of active flavors 
Nf due to heavy-quark thresholds in the b coefficient entering ^^'^(Q 2 ; /XpJ. According to 
this, we set Nf = 3 for /i^ = /i^im whereas for fi\ = /i BLM > //^ in there is no ambiguity by 
virtue of F£' P) {Q 2 ] 

A^blm) — 0- Therefore, the bona fide BLM scale setting reads 

^BLM = max {/i B LM, /Wmin} , (5-14) 

where // m j n will be specified later on in connection with the soft part of the form factor. 



B. ay scheme 

The self-consistency of perturbation theory implies that the difference in the calculation 
to order n of the same physical quantity in two different schemes must be of order n + 1. 



6 Restricting the value of a s does not necessarily limit the quark and gluon virtualities in the Feynman 
diagrams to values for which perturbation theory applies. 
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This means that relations among different physical observables must be independent of the 
renormalization scale and scheme conventions to any fixed order of perturbation theory. In 
Ref. 108| ] it was argued that by applying the BLM scale-fixing procedure to perturbative 



predictions of two observables in, for example, the MS -scheme, and then algebraically elim- 
inating ctjyfg, one can link to each other any perturbatively calculable observables without 
scale and scheme ambiguity. Within this approach, the choice of the BLM scale ensures that 
the resulting "commensurate scale relation" is independent of the choice of the intermediate 
renormalization scheme employed. On these grounds, Brodsky et al. in |45[ have analyzed 
several exclusive hadronic amplitudes in the ay scheme, in which the effective coupling 
av(/i 2 ) is defined by utilizing the heavy-quark potential V(fi 2 ). The ay scheme is a "natu- 
ral", physically motivated scheme, which by definition, automatically incorporates vacuum 
polarization effects due to the fermion-antifermion pairs into the coupling. The fiy scale 
which then appears in the argument of the ay coupling reflects the mean virtuality of the 
exchanged gluons. Furthermore, since ay is an effective running coupling defined by virtue 
of a physical quantity, it must be finite at low momenta, and, therefore, an appropriate 
parameterization of the low-energy region should, in principle, be included. 
The scale-fixed relation between the couplings a^g and ay is given by j4f 



«s(>BLm) = Oty{Hy) 



1 , gvOfy) 8C A 



(5.15a) 



47T 

where 

/4 = e 5/3 /iBLM- (5.15b) 

The scales fiy associated with selected pion DAs are included in Table ITU 

Taking into account Eqs. ()5.15|) . the NLO prediction for the pion form factor, given by 
Eqs. (J4.1|) -([4.8| h gets modified as follows 

(5.16) 

^D,NLO (g2) ^ ^D,NLO (g2)= ^.FG) (g2) + 2 ^LO (g2) _ 

We are not going to present predictions in this scheme using the standard QCD coupling, as 
this would require the introduction of exogenous parameters, like an effective gluon mass, 
that cannot be fixed within the same approach but have to be taken from elsewhere. For 
such an application, we refer the interested reader to the analysis of j45[. The connection 
of [111 to the analytic approach, which we will use below, was discussed in detail in @|. 
Predictions for the pion form factor within the ay scheme will be presented below in the 
context of Analytic Perturbation Theory. 
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VI. STRONG RUNNING COUPLING AND NON-POWER SERIES EXPAN- 
SIONS 



A. One- loop case 



In the one-loop approximation we have a rather simple renormalization-group (RG) equa- 
tion for the running coupling constant: 



d a s {fi 2 ) 
din /i 2 

/3l-loop(/^ 2 ) 



-6n 



47T 



with bo given in Appendix El The solution of this equation has the form 

47T 



a 



(i) 



6 ln(Q 2 /A s 



(6.1) 
(6.2) 

(6.3) 



where A = Aqcd is the QCD scale parameter. A well-known problem here is the appearance 
of an IR pole at Q 2 = A 2 , which spoi l s the analyticity of the QCD running coupling. 

In a series of papers 0, EH llOSl 110l | Shirkov and Solovtsov introduced an analytic 
running coupling that avoi ds by construction the Landau singul arity , thus generalizing earlier 
attempts by Radyushkin [11 if and Krasnikov and Pivovarov [ll2j |. To this end, they used 
the spectral representation for the QCD running coupling a s (Q 2 ) (the bar over a s means 
that the analyticity property is valid) and expressed it in the form 



1 f°° 
a s (Q 2 ) = - da 
^ Jo 



a + Q 2 
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(6.4) 



without subtractions due to the fact that the spectral density p(er) decreases as 1/ In 2 a for 
large a. The corresponding one-loop spectral density reads 



71 



\b J In 2 (a/A 2 ) + tt 2 
and provides the one-loop singularity-free coupling function 



(6.5) 



a 



(i) 



(Q 2 M 2 ) = ^ 



+ 



A 2 



_ln(g 2 /A 2 ) A 2 -g 



(6.6) 



The first term on the RHS expresses the standard UV behavior of the invariant coupling, 
while the second one compensates the ghost pole at Q 2 = A 2 and has a nonperturbative 
origin, being suppressed at Q 2 — > oo. 

Let us now consider powers of the analytic coupling function. By performing an analytic 
continuation of the k-th power of the function ([6.3(1 in the complex Q 2 -plane, one determines 
the corresponding spectral functions p^\a), (k = 1, 2 . . .): 



) Im 

fro 



ln(-a/A 2 ) 



(6.7) 
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which in turn determine the analytic image A^(Q 2 ) of a K s'{Q 



i.e., 



da 



For k = 1, 2, . . ., we have 



-4+i (<2 2 ) = - 

which for k = 1 reduces to 

i(i)/ n 2 



4tt\ ^ ij (Q 



7T 



and ^(Q 2 ) = a£\Q 2 /A 2 ) 



ze 



kb J d\nQ 2 



g 2 A 2 



+ 



ln 2 (g 2 /A 2 ) (A 2 - g 2 



Notice at this point some key properties of these functions: 

• each A^\Q 2 ) with k>2 tends to zero for Q 2 — ■> 0; 

• each A^{Q 2 ) has exactly k — 1 zeros for g 2 G [0, oo); 



when g 2 — > oo, each *4.j^(g 2 ) 



l/ln fc [g 2 ] tends to 0. 



Q 2 ^oo 



(6.8) 



(6.9) 



(6.10) 



These properties are universal in the sense that they do not depend on the loop order. The 
functions ^4fc(g 2 ) are used in the so-called Analytic Perturbation Theory (2I M. la IHR IBtL Eif . 
where standard perturbative series, for example, for the Adler function 



r« s (g 2 )" 


"-J 
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(6.11) 



is recast into a non-power series expansion to obtain 

A(g 2 ) 



D APT (Q 2 " 



1 + 



7T 



7T 



(6.12) 



The one-loop expressions for ,4.1 and 4^ are given in ()6.9|) and ()6.10|) . respectively. 



B. Two- loop case 

In the two-loop case the situation is more complicated. The corresponding /3-function 
reads 

with the first two beta coefficients given in Appendix El Integrating the RG equation (j6.1j) . 
we obtain the transcendental equation 

L Q = -Tp^Tj- - ci In U + -y^V) , Cl = I , I/q = In (g 2 /A 2 ) . (6.14) 
a{L Q )b V a{L Q )b Q J 65 
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As has been shown in |95j, the two-loop running coupling in QCD, being the solution of this 
equation, can be written via the Lambert W-i function 

-l 



a 



(2) 



(Q 2 ) 



An 



b Q d 



(6.15) 



For some more explanations we refer the interested reader to |39j|, Appendix C, Eqs. (C15) 
and (C20) in conjunction with figure 5. By performing the analytic continuation of fu nction 
()6.15j) in the complex Q 2 -plane, the spectral function p^ 2 \a) can be determined |ll3j |: 



P (2 V) 



An 
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where 



z(a 



— exp[— a/ci + z(l/ci - l)n] 



Then, the analytic coupling otPiQ 2 ) in the 2-loop approximation becomes 



da 
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(2), 



cr 



<j + Q 2 
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(6.16) 



(6.17) 



(6.18) 



However, this expression is too complex to be treated exactly. For that reason, Shirkov and 
Solovtsov suggested in [2j to use instead the approximate expression 
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(2,approx) 



(Q 



An 

t(L Q ,c) 



a s [£ {L Q , ci)] 
= L Q + c In 



L% + 47T2 



(6.19) 
(6.20) 



- exp(£) ' 

with the same Lq as in (J6.14j) . This expression reproduces both the UV two- loop asymptotic 
behavior as well as the value at the infrared fixed point Q 2 = rather well. More specifically, 
above about Q 2 > 1 GeV 2 , it resembles the exact result with an accuracy in the range of 99% 
and can be used for all higher Q 2 values. Note in this context that the one-loop expression 



a 



(Q 2 ), Eq. ()6.3|) . can be represented by 



Itt 



afKQ 2 ) = -a s {L Q ) 



(6.21) 



The only feature not yet taken into account in the above approximation is the matching 
at the quark-flavor thresholds: M 4 = 1.3 GeV, M 5 = 4.3 GeV, and M 6 = 170 GeV (with 
Mi = M2 = M3 = 0). However, taking into account this matching, the approximate formula 
f!6.19|) starts to become inaccurate. As a result of the interpolation procedure, we obtain 
then in this (so-called "global" fit in the Shirkov-Solovtsov terminology 0] — abbreviated by 
the self-explaining label "fit") case another approximation: 7 
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(2,fit) 



(Q 2 ) 



An 
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111 ^2 ' C 21 
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v 21 



(6.22) 



7 This interpolation is based upon data contained in 
us by B. A. Magradze. 
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1E3 



and also on unpublished data provided to 
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TABLE III: Parameters entering Eq. 1)6.22(1 for different values of the QCD scale parameter Aq^ d . 



Parameters 


a£d 3 = 350 MeV 


<CD 3 = 400 MeV 


<CD 3 = 4 50 MeV 


M 
A21 

ln^/l GeV 2 ) 


-1.012 
57 MeV 
-5.738 


-1.015 
67 MeV 
-5.412 


-1.091 
69 MeV 
-5.349 



with the parameters and A 21 listed in Table 11111 The quality of this approximation 
ensures a deviation less than 1% in the whole Q 2 interval and is illustrated in Fig. [7| 
To fix the parameter Aq^ 3 , we use |86| 



af\m 2 z ) = 0.120 



that gives us 



A^" 3 = 400 MeV . 



(6.23) 



(6.24) 



Let us now focus our attention to powers of the analytic coupling function. By performing 

the analytic continuation of the k-th power of function (J6.15)) in the complex Q 2 -plane, one 

(2) 

determines the corresponding spectral functions p k (cr), k = 1, 2 . . .: 



Air 
b d 
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1 + Wi(z(*)) 



which in turn provide the analytic images A^\Q 2 ) of a K s\Q 2 ) 
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(6.26) 
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FIG. 7: The solid line shows the approximate expression, given by ai 2 ' fit ^(Q 2 ), Eq. (|6.22|) . whereas the 
bullets represent the exact values of ct~s\Q 2 ) taking into account heavy-quark threshold matching. 



24 



0.001 0.1 10 1000 




0.001 0.1 10 1000 " 5 ~ 4 " 3 " 2 _1 



Q 2 [GeV 2 ] lg 10 [Q 2 /l GeV 2 ] 

FIG. 8: (Left) The solid line represents the approximate expression A^^ (Q 2 ) given by Eq. I|6.28[) . whereas 
the bullets denote the exact values of A^\Q 2 ) taking into account the heavy-quark threshold matching. 
(Right) Comparison of A^^ (Q 2 ) (dashed curve) with [cti ,¥it \Q 2 )] 2 (solid curve). Note the modified scale 
of the abscissa. 



These functions obey a more complicated recurrence relation: (k = 1,2 



dAf\Q 
<91nQ 2 



47T 



I (2) 



and A {2) {Q 2 ) = af\Q 2 /K 2 ). (6.27) 



As a result of the interpolation procedure, we obtain in the "global" case the following 
approximation for k = 2 



4 2 ' fit) (Q 2 ) 



Air 



b (N f = 3) 



1 

L2 



exp(L) 



[l-exp(L)]' 



(6.28) 



with the parameters cf^ and A 2 2 being listed in Table HVl The quality of the approximation is 
high with the deviation restricted to about 1% (10%) for Q 2 > 1 GeV 2 (Q 2 < 0.1 GeV 2 ), as 
illustrated in Fig. [HJ One sees from this figure that for Q 2 > 10 GeV 2 the difference between 
the exact and approximate expression starts to be negligible with the sizable deviation being 
confined in the region Q 2 < 1 GeV 2 . 

C. Factorization of the pion form factor at NLO under analytization 

The analytization procedure of the pion form factor at NLO leads to ambiguities, first 
discussed in j29|. The key question is: according to what analytization prescription are we 



N f =3 



TABLE IV: Parameters entering Eq. Q6.28JI for different values of the QCD scale parameter Aqq D " 



Parameters 


A^- 3 = 350 MeV 


<CD 3 = 400 MeV 


<CD 3 = 4 50 MeV 


c 22 

A 22 

ln(A 2 2 /l GeV 2 ) 


-1.549 
29 MeV 
-7.088 


-1.544 
34.5 MeV 
-6.734 


-1.534 
41 MeV 
-6.399 
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replacing the running strong coupling and its powers by their analytic images? In fact, it is 
possible to impose the analytization of the NLO term of F^ act following two different main 
options: 

• In keeping with our philosophy of the analytization of observables as a whole (Hoi . l6ll | , 

we may adopt a Maximally Analytic prescription and use in the NLO term of the pion 
form factor also the analytic image of a 2 . This amounts to 

^ Fact (Q 2 ;^)] Max An = «i 2) (^)^ LO W 2 ) + ^4 2) (^)^ NLO (Q 2 ;^), (6.29a) 

which will be evaluated with the aid of Eq. (J6.28j) . 

• Another procedure, we call Naive Analytic, replaces the strong coupling and its powers 
by the analytic coupling a s and its powers a 2 everywhere in the NLO term of i^ act . 
This is actually the analytization procedure proposed in and amounts to the 
following requirement 

K act (Q 2 ;^)] NaivAn = -^\tA)^°{Q 2 )+ l - K 2) (^)] 2 ^ NLO (Q 2 ;^). (6.29b) 

Note that the naive analytization does not respect nonlinear relations of the coupling 
owing to different dispersive images. 

Anticipating our detailed numerical analysis of the pion form factor using APT, we define 

AFHQ 2 ) = K aCt (Q 2 )] MaxAn - K act (Q 2 )] NaivAn , (6-30) 
which provides a quantitative measure for the analytization ambiguity. 



VII. PION FORM FACTOR AT NLO: NUMERICAL ANALYSIS AND COMPAR- 
ISON WITH EXPERIMENTAL DATA 

In this section we would like to present our predictions for the pion form factor utilizing 
the BMS pion DA and pQCD at the level of NLO accuracy. First, we consider the standard 
perturbative approach with different scale settings within the MS scheme and continue then 
with a detailed discussion of the pion form factor as a non-power series expansion of the QCD 
analytic coupling. To this end, we employ the analytization procedures discussed before to 
obtain Q 2 F^ act in the MS scheme, with different scale settings, and also in the ay scheme. 
To confront our theoretical predictions with the experimental data in the last subsection, we 
will include the soft non-factorizable contribution, modelled on the basis of local duality. To 
join properly the hard and soft contributions, local duality and the Ward identity at Q 2 = 
will be employed in order to ensure that each of these contributions is evaluated in its own 
region of validity, according to the factorization of the parton and hadron representations. A 
comparison of these predictions with the corresponding ones obtained with the asymptotic 
pion DA will be included. 
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A. Standard perturbative approach 

As outlined in Sec. the NLO prediction for the pion form factor, as any other finite- 
order prediction, contains a theoretical uncertainty stemming from its dependence on the 
renormalization scale and the scheme used. This dependence is, however, reduced in 
comparison with the LO prediction due to the inclusion of the NLO correction. To quantify 
these statements, we plot in Fig. EJa) the ratio R NLO (Q 2 ) = F^ LO (Q 2 ) / (Q 2 ) and in 
Fig. E(b) the result for the factorized form factor at NLO, using the BMS DA in the MS 
scheme with different scale settings. The main observation from these figures is the strong 
sensitivity of -Rnlo(Q 2 ) an d the moderate dependence of F^ act (Q 2 ) on the scale-setting 
procedure adopted — especially at Q 2 values accessible to present experiments. 

Let us discuss these figures in a systematic way. 

• For = Q 2 , the ratio -Rnlo(<3 2 ) is positive, large (on the order of about 50%) and 
decreases very slowly, while a s is small (~ 0.3). As a result, the LO contribution is 
about twice as the NLO one and the form factor is small. 

• Using the FAC scale setting, the whole NLO contribution vanishes, so that also the 
ratio is zero. In this case, the form factor is rather moderate down to momenta of the 
order of 10 GeV 2 , where the QCD effective coupling becomes of order unity. 

• Applying the PMS scale setting, the NLO contribution is negative with -Rnlo(Q 2 ) 
being small and also negative down to a critical value of Q 2 ~ 6 GeV 2 (see Table 
ITT]) , where the absolute value of the NLO contribution becomes equal to the LO one 
and the form factor becomes zero. For this scale setting, already at Q 2 ~ 6 GeV 2 , the 
QCD effective coupling starts "feeling" the Landau singularity and becomes excessively 
large, while above 10 GeV 2 the form factor is rather moderate. 

• Adopting the BLM procedure, the results are quite similar to those obtained with the 
PMS scale setting with respect to the ratio -Rnlo(Q 2 ), whereas the form factor now is 
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FIG. 9: The ratio Rnlo(Q 2 ), (a), and the NLO results for Q 2 F^ act , (b), in the MS scheme with various 
renormalization scale settings. The dashed line corresponds to ^ — Q 2 , the solid line to the BLM scale 
setting with ^ in = /x 2 , = 1 GeV 2 , while the dash-dotted one denotes the result obtained with the BLM (a) 
and FAC (b) scale settings. The analogous result for the PMS scale setting is shown as a dotted line. Note 
that in both panels the BMS DA has been employed and that for the FAC and PMS scale settings Nf has 
been fixed to 3. 
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FIG. 10: The NLO analytization ambiguity Q 2 AF% n (Q 2 ) (left) and the ratio AF% n / [^J act ] MaxAn of the 
NLO analytization ambiguity relative to the factorized pion form factor, computed with the "Maximally 
Analytic" procedure, (right) within the MS scheme with various scale settings: /ip = Q 2 (dashed line), BLM 

(dotted line), BLM (solid line), and the ay-scheme (dash-dotted line). The curves shown correspond to the 
BMS DA. 

negative and very large below 50 GeV 2 (lying outside the range of Fig. Mjo)) because 
the corresponding NLO correction is again negative and even larger. The reason for 
this behavior is that in this scheme the typical parton virtualities in the Feynman 
diagrams are much lower than the external scale Q 2 (see Table |H} giving rise to a 
large value of the QCD effective coupling. 

• The BLM scale setting has two distinct regimes, characterized by the fact that the ratio 
-Rnlo(Q 2 ) changes its sign around 20 GeV 2 : in the regime below this momentum value, 
the result for the form factor resembles that found with the /x^ = Q 2 scale setting, 
though its fall-off with Q 2 is not that steep. On the other hand, above 20 GeV 2 , the 
form factor almost coincides with the one calculated with the PMS scale setting. 

A further complication: it is not clear how to implement quark-mass thresholds when using 
the FAC and PMS scale settings. Therefore, the predictions shown have been obtained by 
fixing Nf = 3. This is because both scales depend on (3q and this induces discontinuities in 
the form factor at the quark-mass thresholds. For that reason, we refrain from using the 
FAC and PMS schemes in our further considerations. To summarize, all scale settings can 
be safely used above about 20 GeV 2 , while at smaller Q 2 values, the PMS and FAC settings 
become unphysical, whereas the BLM and /i| = Q 2 scale-setting procedures can further be 
used at values of Q 2 exhausting the validity domain of pQCD. On the other hand, the BLM 
scale setting remains inapplicable up to scales of the order of 50 GeV 2 (see Fig. EK). As 
already explained before, no predictions in the ay scheme have been shown because this 
would require the introduction of exogenous IR regulators. 

B. Use of non-power series expansions 

We turn now to the results obtained in APT. To exploit the effect of the analytization 
ambiguity on the factorized pion form factor, according to (|6.30|) . we plot in Fig. El (left 
panel) 
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FIG. 11: NLO predictions for Q 2 F^ act vs. Q 2 , using the "Naive Analytic" (a) and "Maximally Analytic" 
(b) procedures and employing the BMS DA. The dashed line corresponds to /i 2 ^ = Q 2 , the dotted one 
denotes the result obtained with the standard BLM scale setting, whereas the solid line shows the result 
calculated with the modified BLM scale setting and the cutoff scale (1q = 1 GeV 2 . The results obtained 
with the ay-scheme are displayed as a dash-dotted line in both panels. 
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and the ratio AF^ n / [^f act ] MaxAn (right panel), employing the BMS DA and the MS scheme 
with different scale settings. Analogous results for the ay-scheme are also included; using 
APT there is no need to introduce external IR regulators. 

Summarizing the results in Fig. the main observations are: 



Fact 

7T 5 



The NLO analytization ambiguity, AF* n (Q 2 ), (left panel) and the ratio, AF^/F^ 
(right panel), the latter being computed with the "Maximally Analytic" procedure 
within the MS scheme with the /x^ = Q 2 (dashed line) and BLM (solid line) scale 
settings, is small and almost scaling with Q 2 above about 10 GeV 2 , albeit in the second 
case there is a sign change around 18 GeV 2 . This is because below this momentum, 
the term F^ 1,FG \ which is negative, prevails, while above that scale the term F^ 1 '^ 
becomes dominant due to — \n(Q 2 / fi 2 ^) — in contrast to the former case in which the 
interplay between these two terms is fixed owing to the absence of the log term. For 
that reason, the "Maximally Analytic" procedure with the BLM scale setting enhances 
the form factor at higher Q 2 relative to the "Naive" one. 

The results with the BLM scale setting (dotted line) resemble those computed with 



-i 2 



the ay scheme (dash-dotted line). In both cases, (/Jr) — a^\fJ^) is large and 

negative (cf. Fig.|H] — right panel), while T^' F is also negative. Hence, the overall sign 
of AF^ n (Q 2 ) is plus because the is absent. However, in the ay scheme the shift 

towards smaller values of the a s argument is much less pronounced and consequently 



the enhancement provided by the use of (/i|J instead of a l f > (/i^) 
(see also Eq. (|5.16|) ). 



is rather weak 



Next, we present the results for the factorized pion form factor derived with APT at the 
NLO level and adopting either the "Naive Analytic" or the "Maximally Analytic" procedure. 
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From Fig. [H\ we see that for both analytization procedures the results for the /i^ = Q 2 
(dashed line) and BLM (solid line) scale settings are very close to each other. Note that 
the ay-scheme yields a similar result (dash-dotted line), but with a much smaller steepness 
of the curve at low Q 2 . On the other hand, the standard BLM scale setting (dotted line) 
produces even an exact cancellation of the NLO and LO terms at the momentum value 
Q 2 « 2 GeV 2 (so too behaves the ratio R N lo(Q 2 ) = F^ LO (Q 2 ) / F^° (Q 2 )— see Fig. 0[a)). 
The origin of this cancellation is, however, purely accidental and unphysical: the BLM 
scale at this point is /Iblm ~ 0-02 GeV 2 with a s ~ 0.75, rendering the pQCD expansion 
unreliable. This deficiency is lifted when applying the 'Maximally Analytic" procedure — 
see Fig. ITTT b). Indeed, such is the impact of the "Maximally Analytic" condition that all 
renormalization-scheme and scale-setting ambiguities are diminished, with all results for 
the form factor almost coinciding, as it is obvious from Fig. ITTT b) . Moreover, from Fig. 
ITTb . we can estimate the effect of varying /i 2 ^ = /i 2 , in the BLM scale setting procedure 
by comparing the BLM (black solid) and the BLM (red dotted) curves. Indeed, //| LM just 
varies from 0.5 GeV 2 (at Q 2 = 50 GeV 2 ) to 0.02 GeV 2 (at Q 2 = 2 GeV 2 ), while the difference 
between these two curves is no more than 10% (using the "Maximally Analytic" condition). 
Actually, for /i 2 ^ varying in the range [1, 0.5] GeV 2 , this difference does not even exceed 
the level of 5%. 

Let us close this discussion with some brief comments on the behavior of the ratio 
-Rnlo(Q 2 )- The message from Fig. [T2] is that, except for the BLM scale setting (already 
discussed), all other scale-settings are not sensitive to the analytization procedure adopted. 
The induced differences are indeed marginal, with -Rnlo(Q 2 ) being positive, large, and prac- 
tically scaling with Q 2 for the /i^ = Q 2 scale setting (dashed line), while this quantity in 
the ay scheme (dash-dotted line) exhibits the same behavior but with the reverse sign and 
having approximately half of its magnitude. The situation for the BLM scale setting is 
somewhat transient between these two options, providing with both analytization proce- 
dures enhancement at the low end of Q 2 and reducing the form factor at Q 2 values higher 
than Q 2 ~ 20 GeV 2 . This effect is due to the (negative) term Fi 1 ' gaining ground against 
F^ that becomes smaller because ln(Q 2 / /i^) is growing. 




FIG. 12: Results obtained for the ratio -Rnlo(Q 2 ) using the "Naive Analytic" (a) and "Maximally Analytic' 
(b) procedures and employing the BMS DA. Notations are the same as in figure ITTI 
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C. Non- factorizable Contribution to the Pion Form Factor 



So far we have discussed only the fa ctorizable part of the pion form factor (cf. ()2.2j) 1 



But as argued originally in Refs. 1181,12! EH US, 

and confirmed latter on in several works, 

for instance, in 21, |25[ 0, 0, 179^ . the dominant contribution at low to moderate values 
of the momentum transfer Q 2 < 10 GeV 2 originates mainly from the soft contribution that 
involves no hard-gluon exchanges and is attributed to the Feynman mechanism. At present 
there is no unique way to calculate this contribution from first principles at the partonic 
level. One has to resort to theoretical models, based on assumptions that attempt to capture 
the characteristic features of nonperturbative QCD. In the present investigation we use the 
Local Duality (LD) approach to calculate the soft contribution, in which it is assumed that 
the pion form factor is dual to the free quark spectral density (1! 24]|, i.e., 



1 fSO fSO 

F^(Q 2 ) = — / / p 3 (s,s',Q 2 )dsds' = 
71 Jtt Jo Jo 

with the 3-point spectral density ps(s, s',t = Q 2 ) given by 



1 + 6 sp/Q 2 



2^3/2 



(7.2) 
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where 



Pais, a', t) 



A (s, s', t) = J(s + s' + i) - Ass' . 



(7.3) 



(7.4) 



Here the duality interval s corresponds to the effective threshold for the higher excited 
states and the "continuum" in the channels with the axial-current quantum numbers. The 
LD prescription for the corresponding correlator ji^] implies the relation 



*o = ^ 2 f 2 . 



(7.5) 



A key issue of the soft contribution is the inclusion of Sudakov-type radiative corrections. 
In |25[ only t he S udakov corrections to the quark-photon vertex were taken into account on 
the basis of |ll7| leading to a reduction of the soft contribution by approximately 6% at 
lo w Q 2 ~ 2 GeV 2 and up to 20% at higher Q 2 values. Just recently, however, it was shown 
in |ll8| that taking into account all radiative corrections to the correlator, the Sudakov 
logarithms cancel out. On the face of this finding, we use in this work Eqs. (|7.2j) - (l7.3l) . 

The soft contribution calculated here is consistent with the result obtained in [29( for the 
asymptotic pion DA on the basis of the soft overlap of the pion wave functions, mode lling 



their kj_ dependence in terms of the Brodsky-Huang-Lepage Gaussian ansatz |119l | and 
using a constituent-like quark mass of m q = 330 MeV. Though the crossover from the soft 
to the hard regime and the asymptotic behavior are strongly model dependent, with the 
mass factor exp (— ^m^jxx) (where (5q is the width of the Gaussian distribution, specific 
for each particular pion DA) playing an important role in tuning this behavior — see |2l| for 
a detailed analysis — the trend at lower values of Q 2 up to about 4 GeV 2 is approximately 
the same. Similar results were also obtained in |l20j ] using a Bethe-Salpeter equation and 
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a constituent-type quark mass of m g = 330 MeV. In both approaches mentioned |2 
the quark mass in the hard part was set equal to zero and the effective QCD coupling 
was assumed to saturate at low Q 2 with a transition scale from soft to hard in the range 
Q 2 ~ 12 - 18 GeV 2 . 
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FIG. 13: Calculation of the soft part of the pion form factor in the Local Duality approach. 
D. Comparison with experimental data 

It is time to step up one level higher and consider the total form factor in order to 
compare our theoretical predictions with the experimental data. So far we have considered 
the factorized hard contribution to the pion form factor only at higher values of Q 2 , where 
pQCD is safe. However, attempting to compute the total pion form factor in the full Q 2 
range, according to Eq. (J2.2j) . we have to combine this contribution with the soft part. Recall 
that we have neglected in the hard-scattering amplitude (i.e., in the parton propagators — cf. 
Eqs. flAl|) . (jA6a|) . and (jA7|) ) all parton transverse momenta against the large scale Q 2 and 
integrated out in the pion wave functions all transverse momenta up to the scale But 
below some momentum scale of this order, these contributions in Th start to be comparable 
(especially in the endpoint region where x — > 1) and, a fortiori, the collinear factorization 
becomes increasingly unreliable. To avoid this happening, we have to restrict the evaluation 
of the hard form-factor contribution to that Q 2 domain compatible with the collinear ap- 
proximation. In technical terms this means that below the scale Sq (the duality threshold) 
we have to switch from the parton representation to the hadron representation according to 
local duality. 8 

As we have seen in Sec. EJ fixing the renormalization scale [i\ in all considered schemes 
entails problems related to the small C} 2 -behavior of the factorizable term of the pion form 
factor: the NLO-term can reach the level of 50% of the LO part, casting doubts on the 
validity of the perturbative expansion per se. In addition, both terms (LO and NLO) 
generate a fast growth of the form factor at small Q 2 , artificially induced by large values of 
the strong coupling and by a l/C} 2 -factor. The origin of this failure, as stated above, can be 
traced to the violation of the collinear factorization approximation, i.e., the resurrection of 
small momenta in Th that have initially been neglected and absorbed into the pion DA. 9 

Hence, it becomes clear that we must correct the factorization results in the low-C} 2 region 




8 A smooth transition from the partonic to the hadronic regime may go via an intermediate constituent- 
quark formation due to QCD dressing. Because there is no unambiguous way to do this, we prefer to 
ignore this regime here (and refer for a discussion of such dressed quarks to 29] ). 

9 Their explicit inclusion would give logarithmic and power-behaved corrections amounting to Sudakov-type 
exponentials containing perturbative |88| and nonperturbative corrections [6o| . 
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in order to ensure that each contribution lies in the corresponding domain of validity. To 
achieve this goal, we need a conceptual framework. 

This is provided by gauge invariance that protects the value of F n (0) by means of the 
Ward identity relating a three-point Feynman diagram at zero-momentum transfer to a 2- 
point diagram. Consider the L-loop approximation in the LD approach. Then, using the 
replacements s — > Sq _1o ° p and ps(s, s', Q 2 ; p^) — > P3~ loop (s, s', Q 2 ; /ifj, Eq. 1)7.2)1 relates 
F^ B (Q 2 ) to the integrated 3-point spectral density p^~ loop (s, s',Q 2 ; p^), which is now de- 
pendent on Recall that the Ward identity links the 2-point (i.e., axial-axial current) 



spectral density p\ loop (s; /ijy to the 3-point (vector-axial-axial current) spectral density 



L— loop . 



s, s', 0; /%) in the following way 



Ps lo ° P (s, s', 0; ^) 



7i5(s 



s')P2 loop (s;/4) 



(7.6) 



Taking into account the LD expression for the pion decay constant, 
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J TV 
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_L-loop 
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P2 _loop (s; pD ds . 



(7.7) 



one finds 



F^(0;p 2 r ) 



1 . 



(7.8) 



The 2-loop approximati on fo r the spectral density, pi ~ loop (s; p^), can be obtained from the 
e + e~ cross-section R(s) |l2l| . because these quantities in massless QCD are proportional to 
each other, so that 



2- loop/ 2 \ 
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1 + ^s(pj) ( 0£ s {pl) 
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^A-^ F + 6o (H_ C(3 ) 



(7.9) 
(7.10) 



where £(3) is the Riemann Zeta function. Then, Eq. 1)7.7)1 yields the following nonlinear 
relation for the 2-loop effective threshold Sq " loop 



2-loop 
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1 + 
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7T 



a s (pi) 
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2-loop 
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/4 



An 2 / 2 (7.11) 



that replaces the standard LD relation, notably, Eq. 1)7.5)1 . Note in this context that the 
effective 2-loop threshold Sq loop should be used only in formulas containing the 2-loop spec- 
tral density p 2 ~ loop (s, s', Q 2 ). Were we in the position to write down the 2-loop spectral 
density p 2 ~ loop (,s, s', Q 2 ) for all Q 2 values, then we would have obtained via Eqs. ()7.2j) - ()7.7j) 
an expression for the pion form factor valid at 0(a 2 ). Instead, we use the leading-order LD 



expression, F^ D (Q ;/%), and add perturbative 0(a s )- and 0(a 2 )-corrections explicitly in 



terms of Fj act (Q 2 ; p 2 ^}. Recalling Eq. ()7.8j) . we then have 

F; act (Q 2 = 0;/4)=0. 



(7.12) 



The next task is to match this low-Q 2 value with the large-Q 2 result of pQCD, Fj act (Q 2 ; p 2 ^). 
The most straightforward way is to adopt Fj act (Q 2 ; p^) at large Q 2 and correct its singular 
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(~ 1/Q 2 ) behavior at small Q 2 by introducing some reasonable mass scale M 1 via the 
replacement 10 

Ff-(Q 2 ;/4) = FAQ 2 ;^)^r - F ^^Mf^- ( 7 - 13 ) 

However, this expression has the wrong limit at Q 2 = 0, so that one needs to correct it in 
order to maintain the Ward identity (WI): 

e ct - WI (Q 2 ; /4) = -FAQ 2 ; /4) HQ 2 /m 2 ) + f.(q 2 ; hDj^^ • (7.u) 

Here the function is some smooth function with $(0) = 1 and z<f>(z) — > when z — > oo, 
introduced to preserve the high Q 2 -asymptotics of Fj act (Q 2 ; The simplest choice for 
$(2) is $(z) = 1/(1 + z) 2 , yielding 

M 2 f M 2 

= (715) 

The scale parameter Mq should be identified with the threshold 2sq" 1o ° p to read Mq = 2sq" 1o ° p 
because Sq loop is the "natural" scale parameter for the 2-point correlator in the pion case, 
while the scale characterizing the 3-point correlator, corresponding to the form factor, is two 
times larger [671 ]. 

In this way, we finally arrive at 

e ct " WI (Q 2 ; /4) = ( 2g2 - lo ? p 2 +Q2 ) 2 ^(Q 2 ; /4) • (7.16) 

We are now in the position to supply an expression for the total pion form factor valid in 
the whole Q 2 range: 

FAQ 2 ; /4) = F^{Q 2 ) + Ff ct - WI (Q 2 ; /4) . (7.17) 

This expression comprises the NLO prediction for the factorized part under the proviso 
of the Ward identity at Q 2 = and the non-factorizable soft part. [Parenthetically, note 
the explicit fj^ dependence of this expression as a consequence of the truncation of the 
perturbative series (see Eqs. (|4.1|1 - ([4.8|) ).] 

Before continuing with the presentation of our final results, let us remark that a similar 
type of matching has been applied by Radyushkin 24] to describe the pion form factor, 
providing the result 

FM) = ^'°'(^) + (a,/,)[l + QV(2. S „)l-' 

1 + a s /7i 



10 One may think of this scale as corresponding to the maximum transverse quark antiquark separation 
bo ~ Mo still accessible to the hard form factor via hard-gluon exchange just before the crossover to the 
nonperturbative dynamics. 



34 



0.6 
. 5 
0.4 
0.3 
0.2 
. 1 



Q 2 F n (Q 2 ) 






j 







1 2 3 4 5 

12 r/~i„w2l 



Q 2 [GeV^ 

FIG. 14: Comparison of the pion form factor calculated by Radyushkin j^J using the LD approach and 
interpolation from low to high Q 2 (dash-dotted line) with our result for the total form factor given by Eq. 
(|7.17() . computed with the asymptotic DA in the MS scheme with the BL M sc a le se tting and using LO APT 
(solid line). The experimental data are taken from (diamonds) and |l22j . |l23| (triangles). 



illustrated in Fig. El (a) (dash-dotted line). In this equation — whic h follows the Brodsky- 
Lepage interpolation formula for the ^-transition form factor |l24j — the first term means 
the soft form factor calculated with the LD approach, while the second one includes the LO 
radiative corrections. It is evident from this figure that Radyushkin's result is very close 
to that given by Eq. (j7.17|) . evaluated with the asymptotic DA in the MS scheme with the 
BLM scale setting using for the sake of comparison LO APT (solid line). 

Employing the above considerations we now present predictions for the total scaled form 
factor vs. Q 2 in different renormalization schemes and perturbation-theory approaches using 
the BMS pion DA. Figure El shows the results for the standard perturbation theory within 
the MS scheme adopting the /i^ = Q 2 (dashed line) and the BLM (solid line) scale settings. 
In Fig. El we present analogous predictions calculated with the APT. In this case, it is 
possible to include results computed with the BLM scale setting and to use the ay scheme. 
We observe from this figure (left panel) that the "Naive Analytization" gives results that bear 
a rather strong scheme and scale-setting dependence. In contrast, applying the "Maximal 
Analytic" procedure, the arbitrariness in the scheme and scale setting is minimized — Fig. 
E^b) being a graphic proof of that. Note that this figure shows also separately the soft part 
of the form factor, displayed in Fig. El and the hard contributions corresponding to the 
various scheme and scale settings discussed above and presented in Fig. ITTT b). 

The phenomenological upshot of our analysis is summarized in the left panel of Fig. IT71 
where we show predictions for the whole BMS "bunch" of pion DAs 0]. The shaded strip 
incorporates the nonperturbative uncertainties related to non local QCD sum rules and also 
the ambiguities induced by the scheme and renormalization scale setting — in correspondence 
to Fig. El Note that the two broken lines mark the region of predictions associated with the 
asymptotic pion DA. These results can be compared with previous theoretical predictions 
and also with further experimental data to be obtained at JLab (see right part of Fig. El 
taken from [3?!). The data points extending to Q 2 of 6 GeV 2 are expectations from projected 
experiments at JLab after the planned upgrade of CEBAF to 12 GeV (we refer to [37| for 
further explanations and related references). 

These striking findings give convincing evidence that the end-point suppressed structure 
of the BMS type pion DA not only provides best agreement with the CLEO and CELLO 
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FIG. 15: Theoretical predictions for Q 2 F 7T (Q 2 ) obtained with the BMS pion DA using standard pQCD 
within the MS scheme and adopting the ^ = Q 2 (da shed li ne) and BLM (solid line) scale settings. The 
experimental data are taken from jig (diamonds) and |l22j |. |l23| (triangles). 

data (cf. Fig. 0}, it also allows to describe the pion form-factor data with at least the same 
quality as with the asymptotic pion DA — as it becomes evident from the LHS of Fig. El 




VIII. SUMMARY AND CONCLUSIONS 

In summary, the key concepts and merits arising from this analysis are as follows. 

• We worked out interpolation formulas for the analytic coupling constant and its ana- 
lytic second power that take into account heavy-flavor thresholds and greatly facilitate 
calculations. This allowed us to develop a theoretical procedure and apply its numer- 
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FIG. 16: Theoretical predictions for Q 2 F 7V (Q 2 ) using analytic perturbation theory and the BMS DA in 
conjunction with the "Naive Analytic" (a) and "Maximally Analytic" (b) analytization procedures. Different 
scale settings within the MS scheme are used: /x^ = Q 2 (dashed line), BLM (dotted line), and BLM 
(solid line). The dash-dotted line represents the prediction obtained with the ay-scheme. Also included 
are the prediction for the soft form- factor part (solid blue line) and below this, the hard contributions 
in correspondence with the predictions for the total form facto r on the upper part of the figure. The 
experimental data are taken from (diamonds) and |l22j . [l23| (triangles). 



36 




FIG. 17: (Left) Predictions for the scaled pion form factor calculated with the BMS bunch (green strip) 
encompassing nonperturbative uncertainties from nonlocal QCD sum rules [lj and renormalization scheme 
and scale ambiguities at the level of the NLO accuracy, as discussed in Fig. El Fhe dashed lines inside the 
strip indicate the corresponding area of predictions obtained with the asymptotic pion DA. Note that this 
strip contains only perturbative scheme and scale ambiguities at the level of the NLO accuracy, calculated in 



APT with the "Maximally Analytic" procedure. The experimental data are taken from [36| (diamonds) and 
[l22| . [l23| (triangles). (Right) Summary of existing and projected experimental data on the electromagnetic 
pion form factor in comparison with the results of various theoretical calculations; figure taken from |37l | 
(see there for explanations). 



ical realization in order to compute the evolution of the pion DA using NLO analytic 
perturbation theory. The hard form factor was corrected at low Q 2 as to fulfill the 
Ward identity and was added to the soft form factor, derived via Local Duality, without 
introducing double counting. 

• On the theoretical front, we found that replacing the QCD effective coupling and its 
powers by their analytic images — a procedure we termed "Maximally Analytic" — not 
only provides IR protection to the coupling, it also helps diminishing the renormal- 
ization scheme and scale-setting dependence of the form-factor predictions already 
at the NLO level, rendering the calculation of still higher-order corrections virtually 
superfluous. 

• From thephenomenological point of view, our most discernible result is that the BMS 
pion DA |l[ (out of a "bunch" of similar doubly-peaked endpoint-suppressed pion DAs) 
yields to predictions for the electromagnetic form factor very close to those obtained 
with the asymptotic pion DA. Hence, concerns that a double-humped pion DA could 
jeopardize the sound application of pQCD are unduly. Conversely, we have shown 
that a small deviation of the prediction for the pion form factor from that obtained 
with the asymptotic pion DA does not necessarily imply that the underlying pion DA 
has to be close to the asymptotic profile. Much more important is the behavior of the 
pion DA in the endpoint region x — > , 1. 

Looking further into the future is yet more exciting. With the planned upgrade of the 
CEBAF experiment to 12 GeV, the pion's electromagnetic form factor can be studied up 
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to Q 2 ~ 6 GeV 2 [371] , providing crucial constraints to verify the various theoretical predic- 
tions discussed here and elsewhere. The apparently good agreement of our results with the 
available experimental data (see Fig. ITTj) is encouraging. 
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APPENDIX A: HARD-SCATTERING AMPLITUDE FOR THE PION'S ELEC- 
TROMAGNETIC FORM FACTOR AT NLO 



In this section we list the NLO results for the hard-scattering amplitude 2 
83l . 84\ 85[, used in our analysis. 

The LO contribution to T H (x, y, Q 2 ; expanded as in (|2.4j) . reads 



where 



2 7tC p 



87T 



1. 



(Al) 



(A2) 



Cf = (N 2 — 1) /2N C = 4/3, Ca = N c = 3 are the color factors of SU(3) C , and the notation 
z = 1 — z has been used. The usual color decomposition of the NLO correction — marked by 
self-explainable labels — is given by 

4 1} (Q 2 ; /4, /4) = C F ^' F) (Q 2 ; /4) + b (Q 2 ; /£) + C G (Q 2 ) , (A3) 
where C G = (Cp — Ca/2) and the first coefficients of the /3 function are 



20 

4C F + — C A ) T&Nf . 



(A4) 



Here, Tr = 1/2 and Nf denotes the number of flavors, whereas the expansion of the j3- 
function in the NLO approximation is given by 



(3(a s ({i 2 )) = -a s (/i 2 ) < 6 



a 5 (/j 2 

4:71 



+ 61 
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(A5) 



With reference to the application of the BLM [62] scale setting in fixing the renormalization 
point, we single out the 6 -P ro P°rtional (i.e., Afy-dependent) term, given by 
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and present the remainder of Th as a color decomposition in the form 



N T 1 

Q 2 xy 



28 

y 



6 - - ) \nx + ( 6 - - * ' ~ ' " ' 

x) v y 



\ny + In 2 (2 y 



O 2 

+21n^(3 + ln(xj/)) 



(A6b) 



(l.G), 



Q 2 



20 

T 



lnx 



In y 



2 In x In y — 2 In x In y 



x y 

+2 In 2 In ]/ + 21nxln?/ — 2(1 — x — y)H(x, y) — 2 R(x, ?/)] 



(A6c) 
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for the color singlet and color non-singlet parts, respectively. For calculational convenience, 
we also supply the sum of these terms (cf. Eq. (|4.7Jl ): 



n(l,FG) 
H 



{x, y, Q\ f4) = C F l£' F) (x, y, Q 2 ; /4) + C G Tg™ {x, y, Q 2 ) 
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where 
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with Li 2 being the di-logarithm function, defined by Li 2 (^) = — J* — , and 
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APPENDIX B: FACTORIZATION SCALE DEPENDENCE IN STANDARD 
PQCD 

Here we e xamin e the dependence of the hard-scattering amplitude (see, for example, 
BQHHIHii). We start with the representation for T H (ji%,fj%), given by (|2~4]l . to get 



t hOf,/%) = <X 8 (fl R ) 



[Bl) 



with T^\x,y,Q 2 ) and T$\x, y, Q 2 ; /4, /4) as in and flU for the LO and NLO, 

respectively. These functions can be represented as follows 



T, 
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(B3) 
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where C\(x, y,Q 2 ) absorbs all other //|- and //^-independent terms from (|A3|) . Using this 
structure with respect to the /if- dependence, we can conclude that 



dT, 



din jip 



T H (j4, /4) ® ^(a s (/4)) - ® ?M/4, /4) + O 



G' 



(B4) 



where V(a s (/i|_)) is the ERBL evolution kernel (|U2|) . Then, the whole derivative of the form 
factor (Q is 

dF i\n£ K) = ^(/ / F)®^(^4)®[n« s (^))-n« s (/iR))]®i(/'F) 
+ 0(« s 3 ) . 

Recalling that in the 2-loop approximation of the standard pQCD 

da s (/j 2 ) \a s (n 2 
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so that 
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V + O(a 2 s ) = 0(a 2 ), (B7) 



:bs) 



Hence, we conclude that at the level of the NLO approximation of the standard pQCD, the 
violation of the factorization-scale independence is one order of a s higher. 11 



APPENDIX C: TWO-LOOP EVOLUTION OF THE PION DISTRIBUTION AM- 
PLITUDE IN STANDARD PQCD 



The pion distribution amplitude (p n (x,p, F ) satisfies an evolution equation of the form 

dip n (x,ii F ) 



din fip 



V(x, it, a s (n F )) ® ip n (u, fi F ) 



(CI) 



where V(x,u,a s (fi F )) is the perturbatively calculable NLO evolution kernel 

OL ct 2 

V(x, u, a s ) = — - V {x, u) + —^Vi(x, u) . 

47T {A7C) Z 



(C2) 



11 Moreover, the dependence on /ip of the NLO prediction for the pion form factor was investigated in 
, where it was found that these results vary only slightly with /ip rendering the "factorization-scheme 
ambiguity" to be small. 
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If the distribution amplitude <p n (x, /ig) is determined at an initial momentum scale /ig (using 
some nonperturbative methods), then the integro-differential evolution equation (jCl|) can 
be i nteg rated using the moment meth od to give (^-(a^/jf) at any momentum scale /ip. The 
one-|l2j| and two- loop jl26l . Il27l Il28| | corrections to the evolution kernel were determined 
in the MS -scheme, but because of the complicated structure of the two-loop corrections, 
only the n umerical evaluation of the (first few) moments of the evolution kernel was possible 
129l 130j ]. However, making use of conformal-symmetry constraints, the complete analytical 



form of the NLO solution of the evolution equation (jOlj) has been obtained [7(1 EH- We 
note that for /ip — ► 00 the solution of Eq. (jCl|) takes the asymptotic form (p n (x, /i| — > 00) = 

Pasix) = Gx(l - X). 

The pion DA can be cast in the form 

tp v (x, f4) = U (x, s; /i F , nl) ® ip v {a, i4) , (C3) 

s 

where the operator U(x, s; /ip, /1q) describes the evolution from the scale /i^ to the scale /ip 
and represents the solution of an evolution equation equivalent to (jClj) . given by 

U(x,s;i4,i4) = V(x,u,a a (j4)) ®U(u,s;/4,l4) ■ ( C4 ) 



dln//p 



It is convenient to express the nonperturbative input DA, (p n (x, /Xq), as an expansion over 
Gegenbauer polynomials C^ 2 (2x — 1) which represent the eigenfunctions of the LO kernel 
V , i.e., 



(p v (x,fj%) = 6x(l - x) 



l + J2' a M)C 3 J 2 (2x 



m=2 



(C5) 



in which denotes the sum over even indices only. The nonperturbative input is now 
contained in the a m (/ig) coefficients. The Gegenbauer polynomials Cn 2 {2x — 1) satisfy the 
orthogonalization condition 



I dx x{l - x) Cl /2 (2x - 1) C z J l 2 {2x — 1) = N n 8 V 
Jo 



(C6) 



with respect to the weight x(l — x), where 



4(2n + 3) V ; 

The moments of the evolution kernel 

M kn [oc s ^l)] = Cl /2 (2x - l)®V(x,y; a 8 (/x|)) V{1 ~ y) C 3 J 2 (2y - 1) , (C8a) 

x y l\ n 

or, equivalently, the anomalous dimensions 

Ikn [0C a Q4)] = ~ 2M kn [<*s{f4)] , (C8b) 

represent the elements of the triangular matrix (k > n). While the LO kernel is diagonal 
with respect to the Gegenbauer polynomials Cn 2 (only the 7^ = 7n°' 1 elements appear), the 
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structure of the NLO and still higher-order kernels leads to the appearance of off-diagonal 
terms in the matrix of the anomalous dimensions (both types of terms ^nn = 7n as well as 
%ni k > n are present). Accordingly, the solution of the evolution equation (jC4j) takes the 
general form 



U(x,s;t4,i4) 

oo 



n=0 

x(l — x) 



Cr(2x-1) + ^ J2 ' 41(^1 ^)C 3 k /2 (2x-l) + 0( 



k=n+2 
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N„ 



Cl /2 (2s - 1) . 



(C9) 



The effect of the diagonal terms 7„ n = 7 n is completely contained in the factor E n (/j%, 
which is given by 



*'da^ M 



The expansion of the anomalous dimensions in terms of a s reads 
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whereas the lowest order anomalous dimensions can be represented in closed form by 
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(n + l)(n + 2)_ 
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with Si(n + 1) = Y^i=i^-h = VK n + 2) + ip(l), while the function ij)(z) is defined as 
ip{z) = dlnT(z)/dz. Since the anomalous dimensions 7„ coincide with the flavor non-singlet 
anomalous dimensions, i.e., the moments o f the split ting kernels in deep inelastic scattering, 
we can use for 7^ the results obtained in 131, 132l |; viz., 
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0, 7 2 



(i) 



830 



N f - 



34450 
243 
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(i) 



31132 
2025 



Nf 



662846 
3375 



(Cllc) 



where Nf denotes the number of active flavors. The nondiagonal matrix elements jkn 
(k > n) manifest themselves in t he d ^ terms of the eigenf unctions expansion and were 
obtained in closed form in 



d kn (A*Fi ^0) - 2 TT Shrill, f4) C kn > 
JVfc 



(C12) 



where 



(0) (0) 
Tk ~ 7n 



Ik in 



2b 



a*Oo) 



2M-l+(7i 0) -7i 0> )/(2bo) 



(C13a) 



12 For Q 2 = 1.7 — 18.5 GeV 2 this number is 4, whereas for still higher Q 2 values, it starts to be equal 5. 
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and 



Cff = (2n + 3) 



- 7 f - 26 + 8CWU„ 2Cp L4 fen - V# + 2) + ^(1)] 



2(fc-n)(fc + n + 3) 



+ 



(n + l)(n + 2) 



(C13b) 



with 



k + n + A 



k — n 



+ 2ijj(k -n) - ip(k + 2) - . (C13c) 



We turn now our attention to the finite-order solutions of the evolution equation (jC4|) . i.e., 
(jCljl . Denoting the formal solution of the LO equation, which contains only the Vq kernel, 
by U LO (x, s; /x f ,/Xq), the corresponding function E n , defined in Eq. (jClOp . becomes 



a s (n\ 



- 7 i 0) /(2feo) 



a 



(C14) 



Analogously, the solution of the NLO equation, containing both kernels Vq and Vi, will be 
represented by £7 NL (x, s; //p, //q). This expression contains contributions coming from both 



the diagonal (E n ) and the nondiagonal (d^n) parts. One finds in the literature two represen- 
tations for the E^ LO (fip, /1q) function. The form which retains the mani fest renormalization- 
group property [E^°(tf, /i|)^ LO (/i|, fi) = E^itf, reads [Hj 



^(/^o) 



cj(n) 



"■v (1) h -v (0) /, 



e n(/4) = «s(/4) 



17i 0) /(26o) 



1 + 



bi 



2&0&1 



u)(n) 



Alternatively, one can recast E^ LO in the form 2^, 13 



(C15) 
(C16) 



-7i 0) /(2&o) 



1 + 



Airbr 



1 - 



2^ 



uj{n) 



(C17) 



which corresponds to a resummation of the leading logarithms associated with the diagonal 
terms, while the subleading ones are expanded with respect to a s . In this work we employ 
(fUToV 3 

Finally, we systematize below some previous results by recasting them in a form that 
is more suitable for practical purposes. As mentioned in previous sections, the coefficients 
a n (/i|) encapsulate nonperturbative information about the binding dynamics inside the pion 
and correspond to matrix elements of local operators according to the OPE, determined at 
som e low -energy scale, characteristic of the nonperturbative dynamics employed [H HH HH 
Il34l Il35l |. To obtain these coefficients at a higher scale, say, /ip, one has to apply LO or 



The expression l|C15|l is obtained by expanding 7„ and (3 to NLO in Eq. (|C10J| and integrating over a s . 
To obtain the form (|C17J) . one expands the integrand in IjClOJI over a s and performs subsequently the 
integration. 
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NLO ERBL evolution. Specifically the coefficients which correspond to the LO evolution 
equation are given by 

«n°(/4) = a°' L °(/4) =a n (fi)El?04,fi), (C18) 
while those corresponding to the NLO evolution equation can be written in the form 

C(^) = < NL °(/4) + ^ar NL °(/4) , (C19a) 

where 

^' NL °(/4) = aM) E^°(^fi 2 ) , (C19b) 

n-2 

«r NLO (/4) = E^^o)^ NLO (^,^)^(^ ; ^) • (ci9c) 

k=0 

We note that using instead of (|C15|) the expression (|C17|) . would introduce only minor 
numerical corrections in the a^ LO coefficients of the order a 2 s (amounting to, for example, a 
1% relative deviation in a^ LO at 10 GeV 2 ). 



APPENDIX D: NLO EVOLUTION BY TAKING INTO ACCOUNT HEAVY- 
QUARK THRESHOLDS 

We describe here the modification of the evolution formulas, presented in Appendix IU1 
due to the inclusion of heavy-flavor thresholds at M 4 = 1.3 GeV and M 5 = 4.3 GeV (with 
Mi = M 2 = M 3 = 0). First of all, for calculational convenience, we limit our study to pion 
DAs that include at the initial scale /i 2 , only two Gegenbauer coefficients (i.e., eigenfunctions) 

<p„(x, nl) = 6x(l - x) 1 + a° 2 Cl ,2 (2x - 1) + a° 4 &J 2 {2x - 1)1 (Dl) 

and rewrite expressions (|G19|) in terms of these coefficients {aS], 04} in the more compact 
form 

4 L °(/^f) = £ 2 (/4,/io) 4 + £> 20 (/4^o) ! (D2a) 
«4 L °G4) = E*(j4, f4) «4 + D42(j4, Mo) 4 + DM, i4) ; (D2b) 

a n>° (/4) = D n4 (f4, i4) Ei(i4i i4) a ° + D n2 (i4, i4) E 2 (f4, iA) a l + d*o(j4> t4) > ( D2c ) 



where 



a s (fi 



Dnk(/4> fi) = »l) . (D2d) 

We start the evolution from the initial scale /Mq = 1 GeV 2 , which corresponds to Nf = 3. 
When //| e [M|,Mf], we need to change our evolution formulas by adopting the value 
Nf = 4. Finally, when /z| > M|, we need to use the value Nf = 5. Besides changing the 
number of active flavors Nf, we need also to match the initial conditions of evolution in 
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each considered region. This generates the following evolution functions (omitting in the 
following expressions the superscript NLO): 



E n (t4,i%;N f = 3), f4<m, 
E n (f4, Mj; N f = A)E n {M 2 , (4) , (4 e (Mi, Mf] , 
EM, Mi; N f = 5)E n (Ml nl) , ^ > Mf . 



(D3a) 



We define in this way the diagonal part of the evolution equation from a fixed initial scale 



(D3b) 



In this way we are able to derive the non-diagonal evolution functions (for the sake of brevity, 
we omit the explicit indication of the corresponding /i| regions); namely, 



£>i0 (Mf, Mo) 



d 42 (i4,i4) 



d 2 ^i 

£>4o(/4 
Ato(/4 



H 2 ;N f = 3), 



Mi;N f 
M 2 ;N f 

l4;N f = 

Mi;N f 

M 2 ;N f 

l4]N f = 

Ml;N f 

M 2 ;N f 



= 4) + E 2 (faM*)D 20 (M*,$), 
= 5) + E 2 (^ F ,Ml)D 20 (Ml,^); 

3), 

= 5) + £ 4 (/4,Mf)D 40 (Mf,^); 

= 4)£ 2 (/4,^) + EM,Mi)D 42 {Ml,^) , 
= 5)S 2 (/4,^o) + EM,Ml)D A2 (Ml,ri) . 



For n > 4 and 



0, 2, 4, we have 



Using these expressions, we can revert to our previous formulas (|C19|) and write 



4) 



a»> NLO (ul) + 



, 2 > . a s(^FJ , ND.XLO, ,2 

4^ — ™ ' 



(D4) 



(D5) 



(D6) 



D nfc (^|,/xg;JV / = 3), 

£ n (/4 , Ml)E-\yil, M|)D nfe (/x|, M|; = 4) + D nfc (M|, //§) , (D7) 
£ n (/4, Mi)E^(ii 2 , Mf ) M|; iV y = 5) + £> nfe (M|, //§) . 



where 



ND,NLO / 



a n {nl) E n {^l, nl) 



n-2 



Ofe(fxg) #fc(/4, /^o)^Sb(/4> ^0) 



(D8a) 

(D8b) 
(D8c) 



fc=0 



with 



47T 

a s (/4) 



D n k{HFi Mo) • 



(D8d) 
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APPENDIX E: EVOLUTION OF THE PION DISTRIBUTION AMPLITUDE IN 
ANALYTIC PERTURBATION THEORY 



The pion DA satisfies an evolution equation of the form 

j: ' 2 F =V(x,u,a s (f4))(g)(f 7T (u,iJ,l), (El) 
ci In /ip u 

with V(x, u, a) having the same functional dependence on a as in (jC2j) . Let us rewrite this 
equation for the coefficien ts of expansion ()C5|) in terms of Gegenbauer polynomials, using 
the notations proposed in |l30l | and linking them to those in |71[ : 



da n (n 



in 



7: 



(0) 



a,(f4) (i) 



An 



7r, 



«n(/4) + n 



4tt 



0<j<n 



2^C (1) M°)-7 (0) 1. 



iV, 



(E2) 



First, we define the diagonal evolution operator E^ h °(Q 2 , /j 2 ,) = e n {Lo) / e n {L ^ with 



dL,, 



57T 



7i 0) + 



47T 



(1) 
i n 



e„ (L 



(E3) 



where L M = ln(jt/ 2 /A|i) and a^ 2 (/j 2 ) is given by ()6.22|) . i.e. 



a 



(2,fit) 



47T 

M3) 



47T 

M3) 



4 2 ' 6t) (^) • 



(E4) 



Then we have 



r-NLO 



exp — 



L 



2& (3) 



26 (3)^ 



[4 2>fit) (^)] 2 ^ • (E5) 



In principle, the lower limit of integration, Lq, can be chosen to be an arbitrary positive 
number, but it is more convenient to set it equal to the average value of L\p F ] under actual 
consideration. Then, we can represent our solution in a modified form — as compared to 
(IUT9|) : namely, 



;NLO 



a n {nl) a s (//|) 



5NLO 



+ 



v ' d- (u 2 u 2 ) aj ^ 

0<j<n C 3 



47T 



(E6) 



The advantage of introducing the same factor e^ LO [L^] for the whole function a^ n (/_i|) is 
that it ensures exact cancellation of the diagonal terms in Eq. (jE2|) . What is left over after 
this cancellation provides an equation for dj !n (fj,p, /i 2 ,): 



d 



a 2 (/4) 



dL 



I'F 



in 



M, Cj 



«.04)v-^ ...2 



47T 



V 'r/- f/; 2 /v 2 W ^ [ — 

/ j u J,mlMF)W Jvl m," -NLO 



j<m<n 



• (E7) 
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In the NLO approximation this expression becomes 



d 



dL 



M. 



26 (3) 



2 e i 



NLO IT 



(E8) 



and its solution is given by 



M 



26 (3) Af' fit) (L 



-NLO / r \ 

. ^ eNLO (L) 



dL . 



(E9) 



Since we do not take into account the nondiagonal part of the evolution equation (see for 
more details in Sec. I1V[ just after Eq. (|4.8jl ). we can use an approximate form of Eq. (jE6|) . 
viz., 



^An^NLO^ = aM yn ^ 



eNLO (L 



(E10) 



where the functions e„ (L MP ) are defined in a two-step numerical procedure: 

1. We determine first by numerical integration of Eq. (|E5|) the functions e^ um (L) for 
L G [0,10] and n = 2,4. 

2. We construct then interpolating functions e^ LO (L) for all functions determined in step 
one. 

In order to obtain the term /i^) in the approximate formula (|4.15|) for the factorized 

form factor, we also need the LO-part of the evolution, namely, 



=LO 



exp 



W AO) 

Jn 



2& (3) 



Af' &t \L)dL 



An: LO/,,2 
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3LO 



CM 



(Ell) 
(E12) 
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